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Chapter 3 . 


INTRODUCTION 


Rising energy prices over the last decade have made 
large wind energy conversion machines an attractive option 
for electric utilities. Still, such machines must be 
carefully designed to compete with the cost of energy from 
other resources. The capital cost, longevity, and 
maintenance costs of the machine determine the cost of wind 
energy produced, and these factors in turn are strongly 
influenced by the dynamic loads and vibrations of the 
structure. Thus, knowledge of the aeroelastic behavior of 
these machines is essential to reducing the cost of the 
energy they produce. 

This thesis examines certain aspects of the aeroelastic 
modeling and behavior of the horizontal axis wind turbine, 
or KAWT. Some problems of HAWT aeroelasticity are simply 
new applications of rotary-wing aeroelasticity, while others 
are uniquely inherent to these systems. Wind turbines 
operate efficiently at relatively large thrust coefficients 
and high inflow angles, and gravity loads are important in 
some analyses. The rotor is subject to a sheared airflow 
due to tlt> earth's boundary layer and the blades must pass 
through the tower wake if the rotor is downwind. 


The first modern generation of large wind turbines is 
typified by the 100 kW NASA MOD-0 series [1]. The 1^D~0 is 
characterized by: 

1} Relatively stiff, cantilevered blades. 

2) Rotor downwind of the tower. 

3} Root pitch change mechanism and highly 
twisted blades. 

4} Yaw drive and brake system to align the 
machine axis with the wind direction. 

This configuration was also characterized initially by large 

dynamic over stresses, but subsequent measures taken reduced 

this problem substantially [2]. 

The 2.5 MW NASA MOD-2 wind turbine [3] embodies various 
advanced features intended among other things to reduce the 
cyclic loading of the blades: 

1} More flexible blades with a teetering hub. 

2) Rotor upwind of tower. 

3) Tip pitch control and less severe blade twist. 

In this machine, the upwind teetering rotor reduces the 
dynamic loads, but the MOD-2 has not been without dynamic 
problems. 

A teetered, tip-controlled rotor has also been fitted 
to a MOD-0 wind turbine for evaluation of these proposals 
[4, 5]. This rotor was tested both upwind and downwind of 
the tower, and further tests with a more flexible tower are 
being conducted. Statistical data on blade bending moments, 
teeter response, and yaw moments are presented in these 



rtfertnces. Furthtr proposals for advanctd HAHTs ineitida 
frst yaw, mora flaxibla towar, and soft driva train to 
isolate tha rotor fron tha ganarator [6]. 

Although HAWT aaroalasticity is a ralativaly naw fiald 
of study, a considarabla body of litaratura axlsts. Briaf 
reviews of the litaratura and methods of HAWT aaroalasticity 
are presented at appropriate junctures in this thesis. A 
good critical review has been presented by FriadsMnn [73. 
See also his survey of rotary-wing aaroalasticity [83* 
Literature in the related field of helicopter aaroalasticity 
has been extensively catalogued in a book by Johnson [93. 
See also his study of proprotor aircraft dynamics [103. 

Wind turbine aeroelasticity is conveniently divided 
into two areas of concern: stability and response* The 
designer's first task is to assure that the machine is free 
of aeroelastic and mechanical instabilities throughout the 
operating envelope, Nonlinear effects are important for the 
stability problem, and the equations of motion must be 
derived consistent with thir fact. Generally, the equations 
of motion are then linearized about some equilibrium state 
of the system. 

The designer's second task is to assure that the 
machine Is structurally sound and will be long lived. The 
response of the wind turbine to the various unsteady inputs 


fmmwL PMH li 

or rooii.^uMiiY 

such S 8 9fsvity forcin9, wind 9ustSi towsr shsdow, and wind 
shear must be calculated. Nonlinear effects are 9 enerally 
less important for the response problem. 

The equations of motion and external load expressions 
of wind turbine systems ace extremely complex. As the 
derivation of these equations proqresses, it becomes 
apparent that some method must be used to sort out and keep 
only the most siqnificant terms. Followin9 the practice of 
rotary-winq aeroelasticity , an assumed ordering scheme for 
the various parameters and coordinates is usually 
established. This ordering scheme is based on typical wind 
turbine parameters and is therefore different from its 
rotary wing counterpart. 

Many aspects of HAWT aeroelastic behavior can be 
studied by modeling an isolated blade with fixed hub. This 
implies that interactions with the tower motion or between 
blades are negligible, an acceptable assumption for the 
MOD -*0 configuration. The isolated blade has important 
degrees of freedom in flapping out of the plane of rotation, 
lagging in the plane of rotation, and pitching about the 
blade axis. 

Much less research has been done in modeling the 
rotor~tower system. The tower may have side~to-side and 
fore^and-Lf t bending, twisting about its vertical axis, yaw 
drive flexibility, and generator drive system degrees of 
freedom. The rotor may have a teeter or gimbal degree of 


freedom. Each blade has flapping and lagging, but blade 
torsion has usually been neglected in examining rotor-toner 
interaction. 

Both hK>0-0 and MOD-2, and most other large wind 
turbines have two blades, a configuration forced by the 
economic considerations alluded to earlier. As a result, 
the equations of motion for the coupled rotor-tower system 
involve periodic coefficients which demand proper 
mathematical treatment. 

There is a temptation to utilize the digital computer 
and numerical methods to construct an all-encompassing model 
of the HAWT system. In this way the issues of nonlinear 
terms, ordering schemes, and periodic coefficients may be 
sidestepped. Indeed, the equations of motion and external 
loads need not be derived explicitly. While this approach 
is useful in the final stages of design, it is unwieldy and 
expensive for initial design or basic research. The 
essential physics of a phenomenon may be overshadowed by 
lesser details of the model, and the source of a phenomenon 
may be untraceable. 

An example of this approach is the MOSTAS computer code 
for wind turbines [11, 12, 13}. This is a very complete 
package which has roots in various helicopter codes, MOSTAS 
has been built up over the years to model most aspects of 
HAWT aeroelasticity and many different machine 
configurations. 
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A general review of MOSTAS has been made by Ougundji 
and Wendell [14]. Based on that study, eight specific 
recommendations were made. Four were suggestions to improve 
the MOSTAS package; four others were rec'^mmendations for 
further basic research. These four are summarized here with 
the original numbers in parentheses: 

1) Develop simpler models to investigate the 
main origins of aeroelastic phenomena (1). 

2) Examine aeroelastic and mechanical 
instabilities more closely, especially for 
the proposed more flexible systems (3). 

3) Study teetering effects and propellor 
whirl type instabilities (4). 

4) Look in detail ^t generator drive train 
interaction with other system components (7). 

The two parts of this thesis address each of these four 

recommendations to a certain degree. The contents are 

summarized below. 

Part I concerns modeling the isolated wind turbine 
blade of a MOD-O type wind turbine. Modeling techniques are 
reviewed briefly, followed by a detailed development of a 
simple three degree of freedom equivalent hinge model of an 
isolated rotor blade for aeroelastic stability. This 
derivation introduces parameters, coordinate systems, and 
concepts of modeling such as linearization of nonlinear 
equations, ordering schemes, and unsteady aerodynamics, all 
of which are common to many of the modeling techniques. A 
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stability st\idy which identifies the important parameters 
and phenomena is presented. Finally, the equivalent hinge 
model is compared to a simple modal model. 

Thus, Part I addresses reconmiendation one, and, to some 
degree, recommendation two cited above. One criticism of 
the computer codes has been their inability to calculate 
stability boundaries [7]. This part of the thesis addresses 
that need with an extensive stability study. 

Part II concerns modeling the coupled rotor-tower 
system of a MOD-2 type wind turbine. Models in the 
literature are reviewed, followed by the development of two 
building blocks for modeling the rotor-tower system. The 
first of these is the derivation of equations of motion and 
external load expressions for a two-bladed, teetering rotor 
on a flexible support. Blade bending modes in flap and lag 
are included, and hub degrees of freedom are used for 
generality. The final equations are in eleven degrees of 
freedom. The second building block is a general harmonic 
balance method to solve systems of second order ordinary 
differential equations with periodic coefficients. 
Stability, steady-state response, and transient response are 
included. 

A simple three degree of freedom model with nacelle 
yaw, nacelle pitch, and rotor teeter demonstrates the use of 
these building blocks. Whirl flutter and divergence are 


Qx«iiin«a, arid the effact of taeter, praconin9, and 
stiffnass is studiad. Transient and forced response are 
calculated for several cases. 

Thus, Part I I addresses mainly recommendations one and 
three cited above. Some aspects of recoi^endation two are 
studied , and tools are developed for recoomendation four. 
Thresher, Dugundji, Hohenemser, and Walton have reviewed the 
state of the art of HAWT structural dynamics analysis tools 
[15]. They reiterate the need for simple models and 
experimental measurements to validate the complex computer 
codes and to foresee dynamics problems. They also cite the 
need to study the effects of teetering and of more flexible 
systems. Dynamics problems encountered in full scale tests 
confirm these recommendations [12, 4, 5; see also 6}. 

In summary, this thesis aims to contribute to two 
general areas of HAWT aeroelasticity : modeMng methodology 

and the understanding of structural dynamics problems. 
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Chapter 2 


MODELING TECHNIQUES 
FOR THE ISOLATED BLADE 


The HAHT blade is essentially a rotating bean which is 
twisted and tapered. Usually, it does not have coincident 
center of gravity, shear center, or tension center, and it 
may have dicontinuities in various properties. In addition, 
the blade may be preconed and its root may be offset from 
the hub. Thus, there are complex structural and inertial 
couplings between torsion, bending in the plane of rotation, 
and bending out of the plane of rotation. 

Koubolt and Brooks derived the linear differential 
equations of such a b&am without precor.e [16]. In 
rotary-wing and HAWT aeroelasticity , it has been recognized 
that some nonlinear effects are important, and several 
researchers have derived related equations which include 
ordered nonlinear terms as well as precone [17, 18, 19, 20]. 
The most important nonlinear effects have been identified, 
although there is a current controversy about 
torsion-stretching coupling [21, 22]. 

The solution method of choice in these studies has been 
the Galerkin modal approach. Hodges and Ormiston developed 
the associated aerodynamic loads for stability of a uniform 
helicopter blade in hover [23]. They solved the equations 
using coupled rotating modes. Wendell developed aerodynamic 


loads appropriata to wind turbine blades and used uncoupled, 
nonrotating modes to examine aeroelastic stability [20]. 
Kottapalli, et al., studied both stability and response of a 
wind turbine blade using uncoupled rotating modes [24]. Of 
special interest is their assertion that the equations for 
stability should be linearized about a time-dependent 
steady-state response of the wind turbine blade rather than 
a time-averaged steady position. 

An alternate, mixed displacement and stress formulation 
of the equations of motion for a helicopter blade in hover 
has been advanced by Stephens, et al. [25]. They solved 
the nonlinear steady-state equations using a collocation 
method, and the linearized stability equations were solved 
by numerical integration techniques. 

The finite element method has also been used to 
formulate the problem. Hodges has developed a method based 
on the Ritz approach [26], and Friedmann has developed a 
method based on the Galerkin approach [27]. Sivaneri and 
Chopra have presented an application of the method to a 
helicopter blade in hover [28]. The nonlinear steady-state 
equations were solved by iteration directly from the finite 
element analysis; the vibration modes were then calculated 
based on this steady-state deflection, and used for the 
stability analysis. Kamoulakos has applied a similar 
technique to HAWT rotor blades [29]. 
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A sinplt nodtl of a rotor blade which nay be propoaed 
is to replace the flexible, cantilevered blade with an 
equivalent articulated blade with springs at the root to 
represent the structural stiffness. This "equivalent hinge” 
nodel has been used extensively to study helicopter rotor 
blades [e.g. 30, 31, 32, 9]. 

Miller, et al., used an equivalent hinge nodel to study 
the aeroelastic stability of a wind turbine blade {33]. 
However, no ordering scheme was used to consistently retain 
nonlinear and higher-order terms. Chopra and Dugundji 
developed a nonlinear equivalent hinge model to study blade 
response and stability [34, 35]. This nodel was 
consistently derived, but ignores center of gravity and 
aerodynamic center offsets, and the feathering axis was 
assumed to be in the plane of rotation. Liebst used a 
derivative of this model to evaluate active control 
strategies for dynamic load alleviation [36]. 

The chapters that follow present a complete derivation 
of the equations of motion and aerodynamic loads, an 
aeroelastic stability study, and a comparison of the 
equivalent hinge and modal models. Nonlinear and 
higher-order terms are consistently derived, cross-sectional 
offsets are included, and the important distinction between 
blade preconing and feathering axis coning is maintained. 
Gravity, wind shear, and tower shadow effects are ignored 
for this simple stability study. 


Chapter 3. EQUATIONS OP MOTION — 

EQUIVALENT HINGE ICDEL 

The equivalent hinge, model is proposed as a simple 
approach to studying the aeroelastic stability of an 
isolated wind turbine blade. The model is useful for 
parameter studies and preliminary design. This chapter 
presents a derivation of the equations of motion without 
aerodynamic loads, which are derived in the next chapter. 

Some nonlinear effects must be retained, but the god 
is to create a simple model and to avoid over dressing the 
simple mechanism proposed. First, the nonlinear equations 
of motion are derived. Second, the equations are linearized 
in perturbations of the blade motion about a steady-state, 
deflected blade position. Finally, the equations of motion 
are simplified by applying an assumed ordering scheme for 
the parameters and the steady-state coordinates. An energy 
approach is followed in the derivation. 

3.1 KINETIC ENERGY OF A GENERAL ROTOR BLADE 

Two coordinate systems are used to describe the 
deformed position of the HAWT blade, as shown by Sketch 3.1. 
The X Y Z system is fixed in space. Tl.e blade is located in 
the X y z system which is rotating at a constant rate 0 so 


12 


er 3DA^ 

YTijAyp 



PAQC « 

QUAunr 


th«t f ■ Ot. Tht blade root la offiat from the hub axis 
and hat a built-in prtcont angle out of the plane of 
rotation. 



Sketch 3.1 Inertial and Blade Coordinates 

The unit vectors associated with the inertial system 
X T Z are X J K respectively, and those associated with the 

A A A 

blade system x y z are i j k. The transformation between 
these unit vectors is 


(3.i) 


^ A ^ AAA 

[I J Kj ■ Cl j 

wh«re [Tj^3 ■ 


» 

eoMp 

0 

•iMp 

cosy 

siny 

0 

0 

1 

0 

-siny 

cosy 

0 

-»‘MP 

0 

CO*^p 

0 

• 

0 

1 


In terms of the blade coordinates and the hub offset, the 
blade position components in the inertial system are then 

X ■ X cos^p cosy - y siny - z sin^p cosy 
+ Sjj cosy 

Y • X cos^p siny ♦ y cosy - z sin^tfp siny (3.2) 
♦ e„ siny 

2 ■ X sin^p z cos^p 

The corresponding absolute velocity components are 

X ■ (x cos/^p " z sin/Jp - Qy) cosy 

- (y -*■ Ox cos^p - Oz sin^p ♦ Oejj) siny 

Y ■ (x cos^p - z sin^p - Oy) siny (3.3) 

- (y ♦ Ox cos^p - Oz sin^p Oe^) cosy 

2 ■ X sini_ ♦ z cosi^ 

P P 

vhere (’) ■ d/dt. 



The kinetic energy of the rotor blade it given by the 
integral over the blade of the kinetic energy of each 
particle 6m 

T - i m(i)^ -► (Y)^ ♦ (2)^3 da (3.4) 

blade 

In terns of the blade coordinates this becomes 

T ■ f { 4l(x)^ ♦ (y)^ ♦ (i)^] ♦ OC(xy “ xy) cos/_ 
blade P 

♦ (yi - yz) sin^ e„^3 ijO^Cx^ cos^yf 

p n p 

♦ y2 ♦ 2 ^ - 2xz sin^ COSyJ 

^ p p p 

+ 2e X cos^ - 2e X sin^ t? 3 ) dm (3.5) 

n pH pH 

This expression is valid for any coordinate scheme used 
to describe the deformation of the blade within the blade 
coordinate frame x y z, for example normal modes and 
generalized coordinates or equivalent hinge rotations. 

3.2 NONLINEAR EQUATIONS OF MOTION 

In the equivalent hinge model shown by Sketch 3.2, 
three degrees of freedom, and three hinge axes, are used to 
represent the HAWT blade: blade pitch $ about the pitch 

control axis; flapping fi roughly perpendicular to the plane 
of rotation; and lagging f roughly in the plane of rotation. 
Furthermore, the rotations are assumed to be in this same 
order. Spring moments about each axis represent the blade 
stiffnesses. ^ 


(side) 


(front) 


Sketch 3.2 The Equivalent Hinge Model 

The deformed blade coordinate system ^ 9 C deflects 
with the blade. The transformation between the deformed 
blade unit vectors 1 y H and the undeformed blade unit 

A A A 

vectors i j k is 

ti j k] • [f r ^31^2^ < 3 . 6 ) 

where [T2] ■ 


cos^ 

sin^ 

0 

cos^ 

0 

sin^ 

1 

0 

0 

-sin^ 

cos^ 

0 

0 

1 

0 

0 

cos^ 

sin^ 

0 

0 

1 

“sin/i 

0 

cos^ 

0 

-sin^ 

cosB 


Finally, the blade position components in terms of the 
deformed blade coordinates are 






X ■ ^cos^cos^ - fcos/sin^ - (sin/ 

y m ^(sin^cos^ - sin/9cos^sin5) 

*»* f(co8^cos0 ^ sin/lsinfsinfi) (3.7) 

- Ccos^sin^ 

z ■ ^(sin^sin^ *»■ sin^cos^cos^) 
if(cos^sin^ - sin>Jsin^cos^) 

- CcosjfsinS 


The associated velocity components are omitted for the sake 
of brevity. 

To simplify the model, the blade is considered to be 
thin and untwisted, that is, approximately in the 4 9 plane. 
Then, C ■ 0 and the following definitions can be made 


j dm * M. 
blade 

/ dm * I 
blade 


j ^ dm » M. e 
bUde ^ ^ 


/ dm • I 
blade r 


f 

blade 

/ 

blade 


1 • Mb«7 

(n dm • 

(3.8) 


Furthermore, it is convenient to lump together the 
blade moments of inertia and products of inertia and make 
the following definitions 


‘ a 



m 

^IIMUTy 


I, • I.cos^^ ♦ I. sin^^ - sin2^ 

1 ^ (F ff 

Ij • lySin^^ l^zon^f ♦ Z^^sin2^ 

1 12 ■ (ly - I^) 8in2^ + l^^cot2^ (3.9) 

©1 • e^cos^ - e^sin^ 

©2 ■ **■ «^COS^ 


Th©s© can b© r©cogniz©d as th© corr©sponding inartial 
prop©rti©s of th© blad© about a sat of axas rotatad an angl© 
f about th© C axis. Not© the useful properties: 


dl ^/df . -21 12 


dl 


/d, - 


21 


12 


dii/d, ■ ii - 


I ♦I » I ♦I 

1 2 ^ y 


d©i/d^ “ "®2 


d©2/d^ - ©1 
(3.10) 


The blade position components (3.7) are substituted into th© 
general kinetic energy expression (3.5), and the above 
definitions are applied to yield 




T • -*• l2>f^ *i(l2 

♦ l^2^osfi fi$ - (Ij^ ♦ l2)sin^ f$ 

0^[Z^^sin/cos9cos^ cos>98in>9p) 

- Ij^sin^cos^p - H^it^9^cos/tsin$] 

♦ OfCdj^ + I^Hcos^cos^cos/p - sin^sin^p) 

* Mj^ej^CjjCos^ ♦ M^e2CjjSin^sin^l 

♦ O^Cl ^(sin^^sin^p - hsin2ficos6cosfi^) 

* I-sin^ - : . ^cosBsindcosfi^ 

2 p X* P 

- Mj^e^ejjSinyJcostf - Hjje2«jjSin5] 

♦ »iO^[I j^(cos^^co8^^p ♦ sin^^cos^^sin^^p 

♦ sin^^sin^^ - 4sin2/ycos^sin2>J ) 

2 2 2 
+ l 2 (cos 9 + sin ^sin 

2 

- I ^ 2 ^sinfisin 26 cos + cosyJsin5sin2>8p) 

- 2Mj^e2eHsin^sin>?p 

+ 2M^e, e„(cosicosi - sin>Jcos^sintf ) ] (3.11) 

b 1 H P P 

As previously stated, spring moments are used to model 
the structural stiffnesses. The strain energy of the 
deflected springs is 

V - * ^9^^ " ^s^^^ (3.12) 

where , k^, and k^ are the flap, lag, and pitch 

equivalent stiffnesses respectively, and the built-in angles 

fi , t , and $ allow for blade "droop”, "trail”, and "pitch 
s s s 
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setting*. To clarify, is the built**in precone angle of 
the pitch axis, and fi. is the built-in droop angle of the 

5 

blade in relation to the pitch axis. 

The equations of motion are obtained by substituting 
the kinetic energy (3.11). and strain energy (3.12) into 
Lagrange’s equations 


J/m . il ^ iV ^ 
dtU?/ fi 


. IT ^ >v * Q 

d.t\hf) if ^f 

3(IX) _ -il ^ iv = 0 

dt\a^/ id id ^d 


(3.13) 


where the Q^s are applied moments about the three axes 
arising from external loads, generally functions of blade 
position and velocity. The equation is 


I + I^2®os^ d - Jsl^sin2^ d"^~ 21^2 2I^cos^ fd 

* 0[2I ^(sin^cos^cos^p cos^sin^p) ♦ 21 ^^^io^cos^p] f 
- Q[I ^(cos^cos/Jp + sin2>JsinyJp - cos2^cos^cos/Jp) 

21 d 

+ ^(sin2^cos^^cos^,Jp - sin2^sin^y5p + cos2^cos^sin2^p) 

+ I “ sin>Jsintfsin2^p) 


+ 2M e ,e (sinicos^ + cosicos^sind )3 k./! 

blH P P ^ 




(3.14) 


The f equation is 
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- (I i' + li2^^ ~ * 2I^eoi/ /♦ 

- n[2I j^Csin^costfcos^p ♦ cos^sin/p) + ^ 

2 • 

* Ol2I ^cos^sintfcos^p * 21^2^^^^ /^sin^p ’*■ >ssin2^os#co8/p} ] i 

* »jO^[ (I 2> <iin48ln2tfco8^^p ♦ cos^sintf8in2/p) 

- 21 ^ sin^^co8^#co8^/p - 8in^tf8in^yfp 

+ is8in2/fco8tfsin2/p)<«’ 2Mj^Cj^ejjSintfsin^p 

♦ 2Mj^e2ejj(cos>Jcos^p - sin^cos^sin^p) 3 •♦■ 

• (3,15) 

The $ equation is 

(I^sin^^J *► I^)^’ *► 1^2^08^ “ (Ij^-'-Ij^sin^ p 

- * 21^2^os^fi f$ 

+ 0[I ^(cos^cos^p sin2^sin^p - cos2^costfcos^p) 

■♦• 21 sin^sintfcos^ ] fi 

X ^ p 

0[2I cos^sin^cos^ * 21. .(cos^^sin^ >ssin2;9cos^co8^ )] f 

* p X fc p p 

■♦■ JjO [- I ^{sin^;Jsin2tfcos^/?p + ijsin2^sin^sin2/9p) 

+ l2sin2^cos^/J I^2^2sin^cos2^cos^>5p * co%fiQ096ftin2fi^) 
-2Mj^ej^ejjSin/Jsintfsin>Jp + 2M^e2e^costfsin^p] k^d 

- k^^ ♦ Q. (3.16) 

Bquatione (3.14*16) are the coaplete nonlinear 
equations for large deflections and vibrations* The 
solution of these equations with all time derivatives set to 
zero is the steady state position of the rotating loaded 
blade, which will be designated fi f $ , 
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3.3 LINEARIZED EQUATIONS 

For investigating small vibrationst the equations of 
motion are generally linearized in perturbations f 0 alMUt 
the steady position /!^ The aeroelastic stability is 
very sensitive to the steady position due to the relatively 
large centrifugal effects [33]. Each equation takes the 
form 


F(^» ft ^t ft fit • Q (3»17) 

Here, F represents the left hand side of the equation of 
motion, and Q the right hand. This form can be expanded in 
linear Taylor's series as follows 

OF/3^')q^ + ... + OF/3>>)q^ +...•»• (3F/3^)^5 

■ Q Q (3.18) 

o 

where the subscript ( denotes evaluation at the 

steady-state position: p » f • f^t 0 • Note that 

the applied moment is expressed as a sum of a steady load Qq 
which is generally a function of the steady-state 
displacements and a perturbation load Q which is generally a 
function of the perturbations. By definition, the 
steady-state equation is 

*o' «o> ■ f‘o- «o> <3-19) 

Therefore the associated perturbation equation is 
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* ... ♦ * ... + (irAi)gi • 8 ( 3 . 20 ) 


Th« perturbation aquations art further sinplified by 
comparing and discarding higher order terms according to an 
ordering scheme which is reasonable for KAUTs* First define 


^0 - 


. k^/o^l^ 




•h 

2 

V. 


•„/L 

(3.21) 


where L is the blade length defined by Sketch 3.2. Then 
assume orders of magnitude 



O 

o 

^S' ^p' ^0' 

OC.**) 


0(.^) 

^0' ^s' ' ®7 

O(i^) 

I* 

0(.2 ) 



These relative orders of magnitude are based on typical HAWT 
parameters. The steady~state deflections depend on the 
other parameters, and their order is determined by examining 
the steady equations (Section 5.1). From another 
perspective, the assumed orders of magnitude define the 
range of validity for the model. 
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The terms in each coefficient of the equations of 
motion are compared to one another, and terms which are 
smaller by one or more powers of t are discarded to simplify 
the equations. Finally, the ordered perturbation equations 
are divided through by 0 X^. In these equations (*) ■ d/dy 
■ (l/0)d/dt. The / equation is 

2[<VV 

• (3.22) 


The f equation is 


* V‘ * * 'V^p>^' 

- Q^/O^Iy (3.23) 


The if equation is 


♦ tT,(v/n) - e 



(3.24) 


These are the important inertial terms 
hinge model of a HAWT blade. 
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Chaptsr 4. AERODYNAMIC LOADS — 

EQUIVALENT HINGE MODEL 

Applied noments about the equivalent hinge axes which 
arise from aerodynamic forces are derived in this chapter, 
based upon unsteady aerodynamics and strip theory* The 
development of linearized aerodynamic loads is similar to 
the develofxnent of equations of motion in Chapter 3, First, 
expressions for the aerodynamic moments are derived which 
are nonlinear in the blade deflections. These expressions 
are then linearized in perturbations about the steady 
position, and simplified by using an assumed ordering 
scheme. Furthermore, the final linearized moments given 
assume quasisteady aerodynamics, a uniform inflow of air, 
ani a uniform, untwisted blade. 

Theodorsen first developed an unsteady aerodynamic 
theory for a pitching, plunging airfoil in i uniform flow 
[37]. His theory used a lift deficiency function to 
represent the integrated influence of the shed wake. A 
rotary wing has a much more 'complex wake structure, but 
Loewy showed that its effect could be included by using a 
modified lift deficiency function [38]. These theories are 
not strictly applicable to a rotor which has chordwise 
motion of the blades, although they have been utilized* 


, , i' ORHINAL PAQK K 

vVi'.'-'- ' qu/ujtv 

^ ^ .•■•' r»‘ k- 

Greenberg developed an unsteady aerodynamic theory for 
a pitching, plunging airfoil in a pulsating flow, thus 
extending Theodorsen’s theory to helicopter rotors in 
forward flight (39]. Hodges and Ormiston adapted 
Greenberg's theory to study the f lap-lag-torsion stability 
of a hovering rotor by using the pulsation to represent the 
lag motion [23]. 

Friedmann and Yuan modified several strip theories for 
use with rotary wings, including Theodorsen's and 
Greenberg's [40]. They compared the various theories with 
and without modification, and studied the quasisteady 
assumption. Recently, Johnson suggested a convenient 
grouping of terms in the lift and pitching moment 
expressions, which is used here [41, see also 42]. 

Pertinent velocity components of the blade axis 
relative to the inflowing air and U^, and the pitch rate 
w , are shown in Sketch 4.1a. Expressions for these are 
derived later. U is the total velocity as shown. Also note 
definitions for the angle of attack «, the blade chord c, 
and the aerodynamic center offset e« from the elastic axis. 
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Sketch 4.1 Aerodynamic Nomenclature 


The noncirculatory lift 4^^^, the circulatory lift 
the drag d, and the pitching moment are defined by 

Sketch 4.1b. In terms of the present nomenclature, these 
distributed loads are 


'no ■ 

- |,.eCtO, - 

where p is the air density, a is the lift curve slope 
yCj^/Ui C is some lift deficiency function, and is the 


— “«■ nut <• 
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profile drag cotfficitnt. The raaultant dlitributtd loads 
of Sketch 4.1c are 


q. • -a 
# ea 

p ■ ^ sin< “ d costf 
f c 

p - --f - Z cos« - d 8in« 

{ nc c 


(4.2) 


All Other components are zero. Since sin« ■ U^/U and 

cos* - U^/U, 

■ i/>ac{ict (ic-e )U, - ic«.U - (Ic^-ice ♦e2)J ] 

7 2^^ 4 4 A 4 4. ^ 32 2 A A ^ 


C^C[-U^ + 

p^ • 5/>ac{CCU^ - (■5C-e^)*^Ujl - DU^U) 
p - ^/»ac{CC-UjU^ ^2®"®A^‘>^;^ 


(4.3) 


- DU^U ^c[-U^ + 

where D « small angles of attack a, U may be 

replaced by whereever it appears here. 


4.1 niLATIVE VELOCITY OF A GENERAL ROTOR BLADE 

The relative velocity of the blade is the difference 
between the blade absolute velocity X Y Z and the inflow 
velccity of air at the turbine disk , which is assumed to 
be in the -Z direction. That is. 


(4.4) 


> 


f • 

“x 


X 


» . < 

Y ► 

“z 
‘ J 
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The absolute velocity components were expressed in terms of 
blade coordinates by equation (3.3). These are substituted 
into equation (4.4) and equation (3.1) is used to transform 
the relative velocities into the blade coordinate system. 

0^ - X - Oy cos^p ♦ 

Uy ■ y ♦ Ox cos^p - Ox sin^p ♦ Oey (4.5) 

Uj - i ♦ Oy sin^p ♦ u^^cos^p 

These relative velocity components, like the kinetic energy 
of a general rotor blade (3.5), are valid for any 
coordinates used to describe blade deformation within the 
x y z coordinate system. 


4.2 EQUIVALENT HINGE AERODYNAMIC WMENTS 

For the equivalent hinge model the position of the 
blade elastic axis is given by equation (3.7) with 9*0 and 
C - 0. 


x ■ ^ cos^ cos^ 

y « ^(sinfcos^ -sin^cos^sin^) (4.6) 

z - ^(sinfsin^ f sin^cos^cos5) 

These are substituted into equation (4.5), and equation 
(3.6) is used to transform the relative velocities into the 
deformed blade coordinate system. The pertinent components 
are 
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^ - ^sin^ 6 ♦ Oft-sin^ sin^ ^ cosi cos^coi5] 

f r r 

* Oe^[cos^cos9 *•* sin^sin^sintf] 

+ Uj. jj[“sin^pCos;Jsin^ + cos^^(cosf sini9 - sin^sin^cos^) ] 

(4.7) 

■ ^cos^ fi ^cos^sin^ B - QeyCos^sintf 
■»• n^[sinyj cos/?sin^ - cos^ (cosfsin^ - sin/ffsin^cos^ ) ] 

V 

+ u^^[“sin^pSin^ + cos^pCos^cos^] 



Sketch 4.2 Resolution of Flap Angular Velocity 

Four angular velocities contribute to the total angular 
velocity of the blade: f is about the C axis; is broken 
down into components about the y and z axes as shown by 
Sketch 4.2; 6 is about the x axis; and O is about the Z 
axis. These angular velocities are transformed to the 
deformed blade axis system using the transformations of 
equations (3.1) and (3.6), and they are added to give 
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The only component required is 

ay ■ Otsin^pCos^cos^ cos/Jp(sin^sin^ + sin/Jcos^cos5) ] 

- sin^ fi + cosA 6 (4.9) 

The work of the external forces can be expressed both 
in terms of the equivalent hinge rotations, and in terms of 
blade motion as follows 


iv 






(4.10) 


where and are the motions of the blade elastic 
axis in the if and C directions, and is the rotation of 
the blade about the elastic axis These two families of 
variations are related in the same way as their 
corresponding velocities. First, set 0*0 and u.„- 0 in 
equations (4.7) and (4.9). 


- ;sin^ e 

ly - ;co5^ * cos/icosf 6 

10 ■ “sin^ ^ * cos/Jcos^ 8 


(4.11) 


Then notice that 
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Sfi • fi 6t 

if » f it 

i6 « i it 


and 


it. ■ 


at 




in ■ 

’ea 


at 

iC ■ 
'ea 

"f 

at 


Substitute these into equation (4.11) and vrite 


^’e=i • ^ " ^sin>5 i$ 

« ;cos^ ofi + ^cos^sin^ i$ (4.12) 

t « -sin^ ij! ■*■ cos^tfcos^ i$ 

These expressions are substituted into the work expression 
(4.10) and coefficients of each variation are separately 
equated to give 


L 

Q- « J C^cos^ - sin^ q ] 6 ( 

L 

^ P. (4.13) 

^ 0 ' 

L 

“ / [^cos^sin^ p - ^sin^ p + cos^cos^ q 1 d^ 

^0 H ^ > 


Finally, the distributed loads (4.3) are substituted into 
equation (4.13). 


- i -j/)acf-jc[^cos^ + (^c-e^)sin^lU^ 

^c((^c-e^)^cos^ **■ (^2®^"|ce^+e^)sin^)ij 
* Cc(^c-e^);cos^ * Ce^)sin^]«^U^ 

- [(C+D)^cos^ - Ce^sin^lU^U^) d^ 


(4.14) 
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- J ^•c{-C(^c-«^)<«^Uj ♦ C<u| - D<a^} d< (4.15) 

• J -/cos^fsin^ * (^-•^)cos^os^]U^ 


♦ ^c((ic«e^)^co8^sln^ - (^c^-ict^'*‘8|)coi^os^]^ 
[C(<|c-e^);cos^s5n^ - (^c^“|cce^^Ce^)coi^cos^ Jt^U^ 

*► [C(^c-e^)^sin4]*^Uy - tCfsinyf3u| ♦ tD^sln>J]uJ 

' [ (C'*^D)/co8^sinf Ce^c08^O8^]lLU^} (4.16) 


Cquation8 (4.14*‘16) are the nonlinear aerodynanic 

moments for large blade deflections (although the angle of 

attack has been assumed small), given in terms of the 

velocities, pitch rate, and their derivatives. With all 

time derivatives set to zero, they are the steady 

aerodynamic moments, which will be designated 

These steady moments are required to solve the steady part 

of equations (3.14*16) for the steady blade position 

/ f 6 , It is of interest that some effects of unsteady 
o o o 

aerodynamics remain in these steady moments due to the 
steady pitch rate OCsin>fp...] in equation (4.9). The 
resulting "apparent camber" . terms are important for a 
preconed rotor. 
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4.3 LINEARIZED MOMENTS 

The aerodynamic moments are expanded in perturbations 
about the steady blade position in the same manner as the 
equations of motion. 

Q * Qq + Q - Qo + (4.17) 

Because of the multiplicity of aerodynamic terms, it is wise 
to establish the ordering scheme for aerodynamic parameters 
now, and apply it to eliminate terms as this expansion of 
aerodynamic moments proceeds. 

The guasisteady assumption will be made, C ■ 1, and the 
blade will be approximated as being uniform along the span, 
with uniform inflow velocity. Actually, only an ideally 
twisted blade would have uniform inflow. A flat blade would 
not, but if 0 is taken as the twist at x/l ■ 0.75, the 

S 

aerodynamic coefficients are practically identical. Now, 
define 

X ■ ^ 

y • pacL^/ly e^ - (4,18) 

Note that here the inflow angle X and the lock number y are 
defined in terms of the hinged length of the blade L. Then, 
assume orders of magnitude as before (Section 3.3) 
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0(.‘‘) 

0<.^) 

0(.‘-‘) 

0(.^) 

ite (4.9), and their 
derivatives are substituted into equations (4.14*16), which 
are then expanded and ordered, and divided by 0 1^ to qive 
the linearized aerodynamic moments 

Q^/0% - 

* - V*o * 2 V' * 

* ^1’’ * 3’’®H ■ 16’'^p ■ * 

- Iy(3x2-8X^ *‘id^)]d (4.19) 

12 ' 00 

V“% ■ ^ 

* tir^J*X-3B^) - i,(lc-5^)(3X-2»^)]r 
* t-|rX(V>«p><3^-2»0> - i-2r<2®-V^-*«o>’^ 

♦ i,(4X-3« )(/» V )? - i^r(«X-3» )« (4.20) 

12 o o p 12' o 

5/0^1, - t-^^„»^(4x-3»_^) - - ir^/ 

. t-lr^,Ux-3»^) - ^,S^(3X-2BJ]f _ 

* '■3-2’’®^ * n’’®*A ■ h*A * 

- i.y(6x2-8X^ +3^2)/- lv(4X-3tf )f 
24^ o o 24' 

* (4X”3^ )]d 

6 A 8"^o 12' o o 




r/8 

X 

z 

®A 


The velocities (4.7), the pitch ra 


(4.21) 


VT' *■ 

These are the important aerodynamic moment terms for the 
equivalent hinge model of a HAWT blade. 




OMOtNW. IMW M 


81 3«)A^ jAHILiPO 
YTI '^Up roo^ 10 


Chapter 5. THE COMPLETE EQUIVALENT HINGE MODEL 


This chapter suounarizes the results of the preceding 
two, and presents the steady-state equations. When the 
linearized aerodynamic moments (4.19-21) are substituted 
into the linearized equations of motion (3.22-24) and 
subtracted from both sides, the complete perturbation 
equations take the form: 





>' 


/ 

/ 

[M] ^ 

r 

i * iC] i 

f 

► Ck3 ^ 

f 


f 

^ i 


t 




{ 0 } 


(5.1) 


The coefficient matrices of equation (5.1) are given on the 
following pages. 
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STEADY EQUATIONS 
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5.1 

The iteady blade poeition is the solution of 

the equations of notion (3.14*16) including the allied 

aerodynamic moments (4.14*16), with all tine derivatives 

zeroed. In their present form the equations of notion are 

quite complex and the steady aerodynamic moments i^re to. 

However, the trigonometric functions can be expanded and the 

ordering scheme applied to keep the most important terms. 

Furthermore, the 6 equation need not be considered. 

o 

This is because 6 ^ will be prescribed for the desired power 

setting, which is related to X and as well as y, c, and 

C. /a. The power coefficient is not constant if any of 
do 

these parameters is varied without adjusting the others. 
Also, y and c are not independent. 


41 


> • MGE I? 

Of Njor QUAirrv 

After some algebra, the ordered equation for is 

found to be quadratic, but uncoupled with 

[2/»p - r8^(8X-3»„)]/»^ 

- * 1 ♦ e^r. - AylBX-SS^j^p - |r(l5-*A)l^o 

* ♦ 

. 0 (5.5) 


The ordered equation for is linear and dependent on 

[>" - (v-»p>'>'o ■ 

'vs * • F^do/* * 

* 2ir(6‘^-8XS-3«2) - jir(V>'p>^‘^''”o>‘ * 


In practice, 6^ is prescribed, and these equations are 
solved sequentially 'rr and After calculating the 
blade steady position, the stability of the blade ia 
examined by extracting the eigenvalues of equation (5.1). 
These occur as complex conjugate pairs, «j t or as real 
roots ay Both the damping a\ and the frequency are 
expressed as ratios per revolution because nondimensional 
time ■ Ot is used. 
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Chapter 6. 


No conplete paraAttric study is attempted due to the 
number of parameters, and the fact that any sp;-cific blade 
design can be easily studied as necessary. Rather, a 
standard case similar to the NASA NOO-0 wind turbine is used 
[1]. Two subcases are studied: a rotor preconed downwind 
and a rotor without precone. Table 6.1 lists the 
parameters. (Tables and Figures appear together at the end 
of this thesis for convenient comparison to one another.) 
The effects of key parameters are examined in relation to 
the standard case by holding all parameters at their 
standard value except the parameters being plotted. 

For the purposes of this section, the blade is assumed 
to be uniform. pitching moment of inertia 
about the center of gravity and ej is the distance the 
cross-section center of gravity is forward of the elastic 
axis, it can be shown that 


■ io * 



3- 

2*1 


where 

T - I /I, - 5, . .,/L 

These definitions allow direct comparision of e^ 
to other studies. 


( 6 . 1 ) 


to e. and 
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Two representations of stability which appear often in 
the literature [e.g. 30, 9] are stability boundaries 

(a . « 0) on the plane of % versus e^ and on the plane of 

3 a A 

versus . Their advantage is that they show the 

interaction of two parameters and much information is 

conveyed by comparing two such plots with one other 

parameter changed. However, a shortcoming is that the 

strength of the instability is not indicated* Such 

information would be needed to assess a faulty blade design. 

Thus, plots of the damping ratio versus the parameter of 

interest, and root locus plots, are both useful as well. 

Some general statements serve to introduce the 
discussion of aeroelastic stability trends and phenosiena. 

1} The flap damping is relatively large. 

2) The lag damping is relatively small. 

3) The pitch damping can approach zero 
under certain circumstances. 

4) Lift couples flap to pitch in a strongly 
unsymmetrir manner. 

The resulting behavior divides the following discussion into 
two parts. First, classical flutter with frequency near v 
or classical divergence may occur. Second, a weak 
instability with frequency near may occur. The flap 

degree of freedom always remains well damped. 
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6.1 CLASSICAL FLUTTER AND DIVERGENCE 

Stability boundaries of the first type on the plane of 
V. versus e_, and on planes of r versus related 

parameters e and ^ , are adequate for the discussion of 
classical pitch-flap flutter and torsional divergence. For 
wind turbines the flexibility of the pitch change mechanism 
reduces the blade torsion natural frequency. The combined 
frequency ratio is important because of its connection 
with the cost of a wind turbine system. 

Figure 6.1 defines the minimum «* required for 
stability of the preconed blade, which increases as the 
center of gravity is moved aft of the elastic axis. Figure 

6.2 is the corresponding plot for the flat rotor. The 

(e^) coupling of pitch and flap in the mass matrix and the 
stiffness matrix can give rise to flutter or divergence 
respectively, although the divergence is not prominent here. 
The motivating force is the large unsymmetric lift term (-|7 
in which couples flap to pitch but not pitch to flap. 

Note that the independent lag instability is insensitive to 
e^ and completely enclosed within the flutter boundary for 
these cases. 

Flutter also occurs when the elastic axis is moved back 

from the aerodynamic center as shown by Figure 6.3. 

Comparing Figures 6.1 and 6.3 shows that for > -.005 or 

i < 0.005 their effects are similar. This again 

A 
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demonstrates that what is important, is the total distance 

that the center of gravity is aft of the aerodynamic center, 

at least for e - e, < c/4, when e* approaches c/6, the 
A I A 

pitch damping approaches zero. This gives the apparent 
asymptote of Figure 6.3 et • 0.0066. However, cases 

such as this do stretch the quasisteady assumption. 

Inspection of the coefficients of the perturbation 
equations (5.1) shows the split personality of which is 
proportional to and moves both the aerodynamic center 

and the center of gravity forward of the pitch axis. Figure 
6.4 shows that, like a fixed wing, sweep forward gives 
divergence while sweep back gives flutter, at least if 
is low enough. 

In all of these cases the flutter boundary encloses the 
boundary of the indepenoent lag instability which proceeds 
at the lag frequency ratio . This suggests the idea of 
separating the three degree of freedom model into several 
two degree of freedom submodels. These would retain the 
steady blade position terms. Their advantage is that the 
new fourth order characteristic equation can be solved by 
hand calculation. The ve'rsus e^ stability boundaries 

of Figure 6.1 were reproduced in Figure 6.5 by this 
technique with the pitch-flap and lag-pitch submodels. As 
could be expected, the flap-pitch flutter boundary compares 
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quite well. The lag**pitch instability boundary is not as 
good an approximation since the system is deprived of the 
flap damping. 
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6.2 LAG XHSTABILXTIES 

Stability boundaries of the second type on the plane of 
versus show the extent of unstable regions and how 
they are expanded by the various parameters, but must be 
supplemented by root locus plots showing the subtle 

interaction of the roots and the severity of the 
instabilities. The pitch frequency is still a key 

parameter, but so are and X and 5^. 

Figure 6.6 clearly shows three regions of instauility 
for the preconed rotor. The familiar pitch-flap flutter and 
divergence is due to the increase in as is reduced 
(see Figure 6.4). .here is a region of flap-lag instability 
associated with the matched stiffness case ■ y which 
may occur when all three frequency ratios are reduced as in 
the case of an overspeeding rotor. The third region, near 
the ordinate, is of most interest. This "stiff 

in-plane" region is very sensitive to the parameters 
mentioned above. 
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The couplings are much more subtle than those of the 
flutter region: 

1) couples f and 6 in the mass matrix. 

2) couples f and 6 in the stiffness matrix. 

3) Coriolis terms couple fi and 

4) There are unsymmetric damping terms due 

to X and 6 . 

o 

5) There is unsymmetric stiffness coupling 

of 4 and 6 due to X and d . 

' o 

The last of these is associated with the torque component of 
the lift which is the prime mover of the lag^pitch 
instability. 

The three degree of freedom character of these 
instabilities is emphasized by another look at the two 
degree of freedom submodels. In Figure 6.7 as before, the 
flutter boundary so calculated compares well with that of 
Figure 6.6. But while the flap-lag region is poorly 
represented, the lag-pitch submodel predicts no boundary at 
all. 

Figure 6.8 corresponds to Figure 6.6 for the flat rotor 
with fi - 0. Especially noteworthy is that the lag-pitch 

Jr 

and flap-lag regions have merged. They are related in 
proceeding at the frequency • 

Figure 6.9 shows the effect of increasing the inflow 
angle X, which greatly enlarges the lag-pitch instability 
region. In this particular plot, the power coefficient is 
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not the same as that of Figure 6.6. A case such as this 
represents a situation where an increased inflow is not yet 
compensated by the pitch setting. This is the only type of 
variation presented in which the power is not held constant. 

Halving the pitch frequency ratio also enlarges 

this region as shown by Figure 6.10. This plot should again 
be compared to Figure 6.6. 

Changing can have a drastic effect on these 

stability boundaries, as Figure 6.11 demonstrates. Here the 
rotor blade has been drooped downwind on the preconed rotor. 
This built-in flap angle has a direct influence on the 
steady flap angle Part of the effects of and of X 

and 6^ also come through the steady flap equation (5.5). 
The leading terms are 

The couplings of f and 6 are all influenced by which 

also increases the Coriolis coupling between f and fi* 

The resulting complicated effect of through is 

shown in Figure 6.12, a plot of the damping of the lag-pitch 

mode a versus the built-in flap angle Had been 

picked as -.15 or -.05, the lag-pitch instability region in 

Figure 6.11 would not have extended to the standard MOO-0 

point at ■ 2.5 and ■ 3.6. The range of which is 

unstable is small for negative fi and the instability is a 

8 

weak one, while positive fi is generally destabilizing. 

o 
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The complex interaction of parameters throu 9 h the steady 
equations thwarts a more specific general statement about 
the lag-pitch instability. 

To examine this instability more closely, root locus 
plots on the iv versus a plane are useful. Only the upper 
half plane will be shown, since the a axis is a line of 
symmetry. The three branches will be labeled conventionally 
as f, and B with respect for their origins, even though 
this may not always represent the nature of the 
corresponding eigenvector. 

The preconed case of Figure 6.12 for the variation 

is replotted in this manner in Figure 6.13. The 
corresponding root locus for the flat rotor is presented in 
Figure 6.14. These plots show the f branch crossing and 
recrossing the a » 0 line in a relatively weak fashion. 
They clearly show the sympathetic participation of the fi 
degree of freedom while the eigenvalue of the f mode is 
dominated by f and 8, 

The migration of the roots as is reduced is plotted 
for both the preconed and flat rotor in Figures 6.15 and 
6.16. In both cases, the 8 branch (flutter) precedes the f 
branch into the right half plane, though only slightly for 
the flat rotor. The 8 branch also finishes near the lag 
frequency while the f branch continues to retreat as v 
is reduced. It is the frequency of this branch which 
coallesces with that of the ^ branch. In fact both the f 
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and $ eigenvalues are dominated by the $ degree of freedom. 
The f instability also can no longer be characterized as 
weak* 

The similar root locus of Figure 6.17 for the preconed 
rotor but with • 2.5 shows the weak nature of the 
flap-lag instability. The stability boundary is just less 
than V ■ 5, confirming Figure 6.10. 
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Chapter 7. COMPARISON OF EQUIVALENT HINGE 

AND MODAL MODELS 

To conclude the discussion of the equivalent hinge 
model, it is illuminating to compare and contrast it with a 
simple modal model derived by this author [20]. For 
convenience, this model is reviewed in Appendix A. The 
equations of motion of an isolated HAWT blade were reduced 
to a three degree of freedom modal model using Galerkin's 
method. One mode each was used for out-of~plane bending 
(w), inplane bending (v), and torsion ( 0 ) , Many details of 
the blade were modeled, including built-in twist, taper, and 
blade cross-section properties all varying along the span. 
All of these entered the squat ions as averages weighted by 
the mode rhape functions. 

The modal model is a three degree of freedom 
mathematical model of the blade. In contrast, the 
equivalent hinge model is a three degree of freedom 
mechanical analog of the blade, and the equations of motion 
derived here are those of this analog. Comparison of the 
two sets of equations reveals many differences, which this 
chapter will discuss. Their sources are: 

1) The deflection shapes. 

2) The coordinates. 

3) Structural couplings. 


52 


ommml nwe n 
W KM qUMJfV 


"^> 4. 1 


Alto, tht modal model it lett general, with a narrower 
ordering tcheme. Xn particular, tquaret and produett of the 
fteady deflect iont tuch at and appear in the 
equivalent hinge equationt, but their eounterpartt do not 
appear in the modal equationt. 

In the modal model, the deflection of the elattic axit 
in the z and y directions, and the tortional deflection, 
respectively, are 

V ■ y^Cx) q^(t) (7.1) 
9 ■ y^(x) q^(t) 

where is the mode shape and q^ is the generalized 
coordinate. Whereas the deflection shapes for the modal 
model were taken as the nonrotating natural mode shapes, for 
the equivalent hinge model the deflection shapes are 


y ■ 

w 




0 < X < L 
elsewhere 

0 < X < L 
elsewhere 


(7.2) 


That it, the blade it a straight line and all of the torsion 
is at the root. Many small differences between these siodels 
arise because the equivalent hinge deflection shapes weight 
spanwise averages of cross-section properties differently. 


I . r 






Modeling the blade with all torsion lumped at the root 
is justified for the common case where the root torsion 
predominates. In other cases a model with the pitch 
equivalent hinge outboard, such as one developed by Chopra 
may be more satisfactory [34]. 

The modal generalized coordinates were expressed as 
sums of steady deflections and perturbations in the same 
manner as the equivalent hinge coordinates. When the above 
deflection shapes were utilized, there is a straightforward 
relationship between the two coordinate sets. For the 
purposes of this chapter, this relationship can be expressed 
as 


and 


Wo 


Vo 




4 ~ 6 
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(7.3) 


(7.4) 


Many apparent discrepancies between the equations dissolve 
when the deflection shapes (7.2) are .Hubstituted into the 
modal equations and these transformations are applied. 

Structural couplings which arise in the modal model 
because of built-in twist and nonuniform cross-section 
stiffness properties do not arise in the equivalent hinge 
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model: Their mechanisms are not present. The most 
important structural coupling is between flap and lag 
bending. The form of the terms in the modal equations 
suggests that an average structural coupling angle can 
be used which is introduced into the equivalent hinge 
stiffness matrix (5.4) as follows 


V/ 

V 


v^cos^d^ * > 1 ♦ ... 

( v^-r?)sin^. costf. * 

^ P boo 

♦ ... 

v^cos^^j^ ysin^^j^ ♦ e 6jj •*■ ... 


(7.5) 


This form was also used in reference [33]. The structural 
coupling angle can be approximated by the blade twist at 
one third span [29]. Similarly, the twist angle at three 
quarters span is generally taken as an approximate blade 
pitch setting 

In principle, this comparision could be carried one 
step fur\;her by using it to relate the integrals of the 
modal model to the parameters of the equivalent hinge model. 
This would perhaps better define these parameters, but if 
such detailed information about the blade is known, it would 
probably be better to use the modal model* 

The two models exhibit good numerical agreement, even 
though overall or typical section parameters are used to 
calculate equivalent hinge parameters. Both models have 
been applied to the MOD-0 wind turbine blade to demonscrate 




this point. The modal model results are taken from 
reference [20]. Figure 7.1 shows the effect of reducing the 
torsional stiffness of the control system, thus reducing r . 
The two models predict the same minimum v required, but 
give divergent results when the pitch frequency a^roaches 
the lag and flap frequencies. Figure 7,2 shows the effect 
of increasing the inflow angle X, again without holding the 
power constant. Figure 7.3 shows the effect of changing the 
precone angle 

The real usefulness of the equivalent hinge model is in 
understanding the effects of the various parameters and in 
simplifying the complex physics of the HAHT blade. This 
model can be used to test concepts, to begin before 
details of a proposed blade are known, or to check the 
results of more complicated analyses. In short, the 
equivalent hinge model is a rotary*'wing counterpart to the 
"typical section" of fixed-wing aeroelasticity . 


Cuapner 6. MODELING TECHNIQUES 

FOR THE ROTOR-TOWER SYSTEM 

The HAWT system is liable to various aeroelastic and 
mechanical instabilities and resonances which involve 
couplings between the main dynamic elements: the flexible 
tower, the yaw drive, the generator drive train, and the 
rotor consisting of several elastic blades and a hub of some 
configuration. Much less research has been presented for 
the rotor-tower system than for the isolated blade. 

Several studies of mechanical instability and the 
effect of static imbalance without aerodynamics have been 
made. Dugundji developed such a model for a two-bladed 
rotor in connection with an experimental study [43]. He 
used an equivalent hinge representation with flap and lag 
for each blade, and two generalized coordinates for tower 
side-to-side and fore-to-aft motion. Sheu used a similar 
but more restricted model with blade lagging and tower 
side-to-side motion only [44]. He investigated ground 
resonance type instabilities for both two and three blades. 

Several studies of aeroelastic stability and response 
have also been presented. Warmbrodt and Friedmann derived 
nonlinear equations of motion and loads for the rotor-tower 
system [45]. Galerkin’s method was applied to study a MOD-0 
type wind turbine with two blades. They used one lag and 


one flap mode per blade, and one mode each for tower 
torsion, side-to-side bending, and fore-to-aft bending. 
Hultgren and Dugundji developed a similar linear model for 
the three’bladed case [46]. They studied mechanical 
stability and forced vibrations due to imbalance, gravity, 
and wind shear. Thresher, et al., examined the response of 
a three-bladed rigid rotor on a flexible tower to 
atmospheric turbulence [47]. Bousmann and Hodges presented 
an excellent experimental study of the aeromechanical 
stability of a three bladed hingeless rotor on a flexible 
pylon [48]. Bundas and Dugundji have conducted some 
experiments on the yaw behaviour : a model wind turbine 
with two blades [49]. 

Recently, Janetzke and Kaza presented an analytic 
rotor-tower model with a two-bladed teetering rotor 
applicable to the MOD-2 [50]. They used a teetering rotor, 
one flap mode for each blade, and a kind of equivalent pivot 
model of the tower nacelle with yaw and pitch degrees of 
freedom. Whirl flutter was investigated by numerically 
integrating the equations in time. This is quite similar to 
an approach ust.i by Hall to study whirl flutter of a 
teetering proprotor [51]. 

Finally, it should be noted that various computer codes 
have been applied to the rotor-tower problem [e.g. 2, 12]. 
However, documentation of the theory used is generally poor, 
and very few parameter variations are given. 
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In the chapters that follow, a linear aeroelastic 
modeling methodology is developed for the MOD-2 type wind 
turbine. Equations of motion and loads are derived in 
closed form for a two-bladed teetering rotor with elastic 
blades on a flexible support. One lag and one flap mode are 
used for each blade, and six general hub degrees of freedom 
are assumed. A solution method for the resulting periodic 
coefficient equations is presented which is applicable to 
stability, steady-state response, and transient response 
calculations. 

The methodology developed is demonstrated with a simple 
yaw, pitch and teeter model similar to that of Kaza and 
Janetzke. A limited study of the effect of imbalance is 
made. Whirl flutter and divergence, as well as other 
instabilities are examined, and the effect of teeter, 
precone, and support stiffness are discussed. Some 
steady-state and transient response results are presented. 

Thus, while Part I mostly concerns the MOD-0 type wind 
turbine with cantilever blades. Part II mostly concerns the 
MOD-2 type with a teetering rotor. 'lutter, divergence, and 
lag instabilities li’e those discussed in Part I are 
possible for the teetered rotor as well, but are modified by 
the interaction of the blades. For the helicopter case, 
Shamie and Friedmann have analyzed the flap-lag-torsion 
stability of a teetering rotor on a rigid support [52]. 
This problem is not addressed by this thesis. 
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EQUATIONS OF MOTION — 
ROTOR-TOWER MODEL 


A mathematical model of a teetering rotor on a flexible 
support, once derived, is a valuable building block for 
investigating HAWT aeroelast icity . The develofxnent 
presented here allows flap and lag modes for each blade, 
teetering motion of the rotor, and six general hub degrees 
of freedom — eleven degrees of freedom in all. Thus, the 
new information given by this derivation is essentially a 
description of the interartion between rotor modes and hub 
motion. The hub degrees of freedom can then be used to 
match the rotor to any kind of tower or support model, from 
simple to complex. This approach is taken by other analyses 
as well [10, 11]. 

The equations of motion are derived in this chapter 
with gravity loads but without aerodynamic loads. The 
latter are developed in the following chapter. An energy 
approach similar to that of Part I is followed, but only 
linear terms will be obtained for the rotor-tower model. 

9.1 COORDINATE SYSTEMS 

Three coordinate systems are used to describe the 
deflected position of a rotor blade as shown by Sketches 9.1 
and 9.2. The inertial coordinate system X 7 Z is fixed in 
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A A A 

space and has unit vectors I J K. The shaft coordinate 

, A ^ A 

system 7^ (Xg Jg Kg) defines the deflected position 
of the shaft axis and hub, but does not rotate. Finally, 

A A JL 

the blade coordinate system x y z (i j k) locates the 
rotating, teetering blade with x as the unbent position of 
the elastic axis. Elastic motion of the blade is described 
within this blade coordinate system. 



Sketch 9.1 


Inertial and Shaft Coordinates 


1 
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Originally* the hub is at the Z 7 Z origin* and the 

shaft axis is at the Z axis. The displaced hub is located 

by three Cartesian deflections g^* q^* and q^, as shown by 

Sketch 9.1. The shaft axis Z is deflected in three 

s 

rotations ^^* ^„* and about Z* 7* and Z* respectively. 

XX z 

Also, the rotor spins about Z^ at a constant rate 0. This 
is not a restriction on the model however, since 
perturbations in rotation speed can be included in 


X / 


Blade 1 




Blade 2 — ^ 


(side) 
Sketch 9.2 


(front) 

Shaft and Blade Coordinates 
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Blade one is arbitrarily chosen as the reference blade 
in Sketch 9.2. The azimuth angle of blade one is y ■ Ot. 
It should be oted that equations of motion and loads need 
be derived only for blade one. The results apply to blade 
two as well, but with |^ ■ Ot ^ » and ^^2 ■ “^tl 
contributions from both blades are then summed as a final 
step in the derivation. The subscripts 1 and 2 are dropped 
from this point except where required for clarity. 

The blades have an instantaneous teeter angle in 

opposite directions, and a built-in precone angle in the 
same direction as shown by Sketch 9.2. Both and are 
positive in the upwind direction for blade one. 

As previously mentioned, elastic motion of the blade 
occurs in the x y z coordinate system: Flap bending w is in 

the z direction and lag bending v is in the y direction. 
The blade is also foreshortened by the bending, which gives 
rise to a deflection u in the x direction (not shown for 
clarity). Torsion and blade stretching are ignored in this 
analysis, because they are generally modes of much higher 
frequency for HAWTs. 

To simplify Sketch 9.2, a- small offset was omitted. A 
preconed rotor as shown is not balanced in teeter and flops 
forward under its own weight. Practical preconed rotors 
might counteract this by using a small undersling e as shown 
in Sketch 9.3. Thus, the teetering axis is actually 

situated in blade coordinates at (e sin>9 , 0, e eosfi ). 
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Sketch 9.3 


Rotor Undersling 


To sununarize, the displacements and rotations in going 
from the inertial system to the blade system arer in order 
of appearance: 

Three Cartesian deflections Qy 
Three rotations 
Shaft rotation f. 

Precone and teetering 
Undersling e. 

Elastic blade deflections u v w. 

The rotations are taken in the order ^y, then then 

and ft then and Thus, the transformation 

between blade unit vectors and inertial unit vectors is 
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[i j k] • [I J K][T^] (9.1) 

where [Tj • (Tyl (T^HT^HTg] • 


— 

» < 

» • 

|. • 

• 

cos^y 0 sin^Y 

10 0 

COS'F -sin'F 0 

cosB 0 -sinB 

0 10 

0 COS^Y 

sin'F cost 0 

0 10 

-sin^ 0 cos^Y 

» * 

0 

r-t 

O 

O 

sinB 0 cosB 
» ■ 


here^ 'Y • f * ® convenience. 


9.2 KINETIC ENERGY 

With all this information, the radius vector in the 
inertial coordinate system can be written for a point x on 
the rotating, deflected blade. For blade one, it is 


A * ^ . M 

[I J K] <Y 


Z 
^ J 


{ \ 

Q, 


n ^ « I 

[I J K] < q 


‘ZJ 


I ♦ (i j ic]^ 


u + X - e sin^ 


w - e cos^ 


P J 


And, by substituting equation (9.1) in, this becomes 

fu X “ e sin^J 


A A A . 
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r ^ 
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► * [T ]i 
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1, 

V 



w - e cosfi 


(9.2) 


(9.3) 
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The corresponding velocity vector is 
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V - [I J K> 

1 

’x 


u^x-esini 
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u 

# 

V 
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w - ecosi 

'"pj 


• 

w 



Recall that blade torsion is neglected in this 
analysis. Torsional moments are also assumed to have little 
effect on the blade bending or the other degrees of freedom, 
particularly when contributions from both blades are added. 
Specifically, torsional moments which arise from inertial 
sources are neglected by assuming that the blade is a long, 
slender beam with all its mass m(x) concentrated along the 
elastic axis x. The kinetic energy of the blade is then 




HJ mV*V dx 
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(9.5) 


Equation (9.5) must be expanded by substituting the 
coordinate transformation (9.1) into it. Kinetic energy 
terms which are quadratic in the displacements produce terms 
which are linear in the final equations of motion. Elements 
in the transformation matrix, *its time derivative, and their 
products must be kept to adequate order to retain all 
quadratic terms, but need not be kept beyond that. These 
matrices are relegated to Appendix B. 

A simple modal model is assumed for inextensional 
bending of the blade, as follows 
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w ■ y^(x)q,,(t) 

V • (9.6) 

L. 2 , , 2 

u - -*««v<Iv ■'>*W<J« 


where 


X ^ X i 

«V ■ ^ ^Tv> " 0 

This form of elastic deformation and the various 
matrices are substituted into equation (9.5). With all 
quadratic terms, the kinetic energy of blade one is 

- S*[qjj(/j^coS)> - j»y) + ^^(/^sinf - fy^)] 

+ S^cos^pt^Jj. + fy^siTif - f\cosf]qj 

+ OS*CqjjSinj^ - q^cosf]/i^ - OSj^cos^p (q^^siny - q^cosf) 

-► nSjjCOS^p[^j^cosf + ^^sinflq^ 

- OS^IqjjCOSf + q^sinfjq^ •»■ OS^sin^pCq^^siny - q^cosfjq^ 

- S^[sin^p(qj^cosf ♦ q^sinf) - Q^cos^plq^ 

2/j^^‘jjSiny - 2^‘^^\cosr ) 

* 2fijy^%inf - 2^^;^cosf ♦(j»j^)^sln^f ■*• 

- 2j»jjy*„sinfcosf )^cos^yl 

- S^CqjjSinj^ - iyCOSK^q^ + QI**(/j^coSf + /^siny]^^ 

- OIj^cos^,tfp[(j?jjCOSy + f^%iTif)fi^ * f^fy) 

■*• OI*^(^'j. ■*■ fy9iT\f - ^'^cosflq^ 

- OI*co5^p[jJjjCOSf ♦ /^sinf] ♦ W^I^cos^^pCl - ^^3 + .•• 
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♦ Ol*^sin>Jpt^*jjCosr ♦ ^\»inr]q^ ♦ ^ 

- OI^^coi^^ptjTjjCOSf ♦ #\iinflq„ “ O^ljcof^p 

♦ ♦ ijO*l^iin*^l, q* ♦ 0*l‘^ftn^j, 

- 0*l_^^iin/^eoi/^ - I'^J^'^eoif ♦ ^\»in|>]q^ 

♦ ♦ ^‘j(»ln|» - ^’^coiflq^ - 

+ OI„eoM^ q^ - 01*^ ^^q, ♦ - i*!-’ 

- OSj^eos/ [q^cosp ♦ q^iinplpj - 20I*eot4|, 

- S^co8>?p[qjjSinf - q^cosj^J^^ "" **^b®®* 

+ 0l*co8^p[^‘j^sinf - f^co9f\f^ * 0I^co8^/p 

- I*co8/pC^’x®®®>' ^ * ^xv^®*^p 

- 2ni 9infi cosfi iq 

xw p p z^w 


( 9 . 7 ) 


In thi8 equation, the following definition8 are uaed: 
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These latter definitions reduce to zero for the casl^'* with 
neither precone nor undersling. 


9.3 GRAVITY POTENTIAL ENERGY 

In this analysis gravity is assxuned to act in the 
negative X direction. The effect of built-in rotor axis 
tilt is not included. The potential energy of the blade is 

L 

V - Xmdx (9.10) 

Prom equation (9.2) and equation (9.1), the height of a 
point X on the blade is given as 

X ■ Qv '*■ tcos^ coeBcos’P sinf sinf coaBsinSr 

A X X A 

+ cos^ sinB)(u+x-esin^ ) 

X p 

+ t-cos^Y^int sin^Y*i"Px^'* 

+ t-cos^ySinBcos'P - sin^^sin^^sinBsin'P 

sin^^cosa^cosB] (w-esin^J ) (9.11) 

I X p 

As with the kinetic energy, the potential energy is expanded 
in terms up to quadratic in the displacements and 
velocities. Thus, the potential energy of blade one is 

V - gM^q + gs^eoi^plco*, - #^sinf cos, + 

-*i9c08/^eo»,[H^qJ ♦ * 9S^t»in, ♦ ,^co»fJq^ 

♦ gS^[cos^p(,»^ ♦ ^^cos ) - sin^p(cosf - 

+ gS*[^Y “ (9.12) 
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The spanwise integrals in this equation are defined as 
before (9.8-9). 


9.4 STRAIN ENERGY 

The elastic model <ji the blade is taken directly from 
reference [20] and is reviewed in Appendix A. However, only 
linear lag bending and flap bending are used; torsion and 
all nonlinear terms are ignored. A teeter spring is 
included, ..nd structural coupling due to built-in twist of 
the blade (x) is allowed. Sketch 9.4 shows the principal 

b 

axes defined by 

b 



V 

Sketch 9.4 Built-In Twist and Principal Axe< 

The strain energy of the blade- is 

0 • ♦ w,/, (9.13) 

Here, is one half of the teeter spring rate, and the 
stiffness coefficients are defined as follows 
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S (EI^ - Zl ^)sin0^cosB^ (yjyj)dx (9.14) 


where 


5 (El2Sin^d^ •*■ El (y^)^dx 


El 


J E,‘ 


d,dC 


and 


x-sect 


El, - f EC^dfdC 


x-sect 


9.5 EQUATIONS OF MOTION 

With the kinetic energy (9.7), gravity potential energy 
(9.12), and internal strain energy (9.13) of the blade in 
hand, Lagrange's Equations are used to develop the equations 
of motion 


_d/3T \ il . ^ ^ 

dt\3^; “ aq aq ~ aq 

n n n n 


®n 


(9.14) 


where the q^ are the generalized coordinates q^^, , . . . r 

q and the Q are corresponding generalized loads, 
w n 

Lagrange's equations are applied first to blade one 
with q q , and q * q in T, V and U. Then they are 

V vl w wl 

applied to blade two with ■ q^^, q^ ■ 

f * f "*■ n , The contributions from both blades are added to 
give the equations of motion in eleven degrees of freedom: 


* OUMunr 


six hub deflections 

and 

rotations. 

teeter, lag bending 

of 

each blade, 

and 

flap 

bending 

of each blade. 

No 

contributions 

from 

the tower or other system components 

are 


yet included. 

Some sinplif ication is wrought by utilizing elastic 
nodes for the conplete rotor instead of separate nodes for 
each blade. Consider the transformation 


s "l ^2 

* '’wj 

a 1 2 


(9.15) 


Then g indicates motion symmetric about the rotor 

WS 

centerline and indicates antisymmetric. Opposite 

definitions must be used to give lag modes which are 
symmetric CC" shaped) and antisymmetric ("S” shaped) in the 
usual sense. 


g ■ ^Cg * a ] 

^a '^l ^2 

g » JjCg - o ] 

''s ^1 ^2 


(9.16) 


These sets of eguations 'are multiblade coordinate 
transformations of a kind for two-bladed rotors. Their 
implications are made clear in Chapter 11. 
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The vector of generalized coordinates and the vector of 
corresponding generalized loads are defined as 


'•’x' 
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X 

<Jy 
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Pz 



Ox 
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«Z \ 
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^a 
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^s 
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q 
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p 
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The equations of motion in these variables take the standard 
form 

[M]{ql + [CHq} + [KHq} - {Q) (9.18) 

The periodic~coef f iclent matrices of these equations are 
given on the following pages. For convenience in comparing 
them to the aerodynamic derivatives of Chapter 10, they are 
also given in individual coefficient form. These equations 
are left dimensional until the tower contributions can be 
added. 
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Chapter 10. AEIK)DYHAMXC LOADS - 

ROTOR-TOWER SYSTEM 


Applied forces and SK»ents which arise froa aerodynanic 
forces on the blades are derived in this chapter in parallel 
fashion to Chapter 4. As with the rotor-tower equations of 
motion, torsion and torsional moments are neglected. This 
assumes that the aerodynamic center is at the elastic axis. 
The quasisteady assumption is made and apparent mass effects 
are ignored. However, preconing, blade twist and taper, 
wind shear, and crossflow over the rotor are included here. 
As in Chapter 9, the contribution of one blade is derived 
first. Then, loads from both blades are summed and the 
symmetric and antisymmetric coordinates are used. 

A cross-section of the deflected blade is shown in 
Sketch 10.1a which defines the deformed blade coordinate 
system f f C. The corresponding unit vectors are f y «. 
The blade has a built-in pretwist d. (x), followed by 

D 

deflections v and w as before. 
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Sketch 10.1 Aerodynamic Nomenclature 

Pertinent velocity components of the blade axis 
relative to the inflowing air and , and the pitch rate 
are reviewed in Sketch 10. lb* The distributed loads of 
interest here, p and p as shown in Sketch 10.1c, are 

7 $ 

distilled from equations (4.3) with the simplifications 
outlined above. They are 


p - %/»ac{CC-U^Uj (10.1) 

where f is the air density, a is the lift curve slope, c is 
the blade chord at station x, C is soae lift deficiency 
function, and D ■ C^^a. 

It is actually more convenient to vork with velocity 
and distributed load coo^nents in the x y z system for this 
development. These are shown in Sketch 10.2a and b, 
respectively. Note that the pitch rate is retained in 
the deformed system; its transformation produces no 
simplification. 



Sketch 10.2 Components in the x y z System 


The transformation between the deformed blade unit 
vectors i y » and the undeformed blade unit vectors i j k is 
taken from reference [20]. Xt is 


where [T.] 

D 


m 

f 4 i 4 ^ 

l« y «J ■ 
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4/4 
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- w'sin^. 
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cos$. 
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-sind 

v’sind^ 

b 

+ w*cos 5 ^ 

b 

sin^. 

b 

cos 5 


( 10 . 2 ) 


When this transformation is used for the velocities and 
distributed loads in equations (10.1), the distributed loads 
become 


[-(C+D) sin^jjCOS^jj w’ - (Csin^^jj-Dsin^^i^ )v* luj 
♦ {C-^D)[cos2^. w' + Sin2^. v‘]U U 

D D y z 

* [ (C+D)sind. cos^. w' - (Ccos^^. -Dsin^tf. )v* ]oJ 

D D D D 2 

- i 5 cC[sintf. v’ cos^. w']i»,U 
b b f y 

JscCCcostf. v’ - sin^. w']» U ) (10.3) 

b 0 ^ z 

p ■ ij^ac{ 2[ (Csin^tf. ‘►D)v’ - Dsin^.costf. w’jDD 
y ^ b b b X y 

- 2[Csin0^costf ^ v' (C+Dsin^^. ) w' )U U 

b b b X z 

- Dcostf, - (C+D)si:n^, U U ♦ Ccostf^ 0^ 
by b y z b z 

♦ jjcCtsin^ v' ♦ cosd w’]« U -ijcC *^U } (10.4) 

b b ^ X r ^ 
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• »s/k:c{ 2t-Cfin^^coi^jj V ♦ (C^Dcoi*^,,>w* ]O^Uy 
♦ 2[(Ccos^^jj+D)v* ♦ Dfin^j^cof^^ 

♦ C%in$^ Uy - iC*D)cot$^ OyO, - Dllntf^ 0^ 

- >*eC[co8#^ V ♦ iln^^ '**J*^®x '*’ **cC*|®y J (10.5) 


Thus, th« velocity components U^, end the pitch 
rete ere required to tormulate the distributed loads. 
For this purpose terms need be kept only up to linear in the 
displacements. 


10.1 RELATIVE VELOCITY OF THE BLADE 

The relative velocity of the blade U is the difference 

between the absolute velocity of the blade V and the 

velocity of the air at the turbine disk. An inflow velocity 

u, in the “2 direction and a crossflow velocity u^^ in 
in cr 

the 'T direction are assumed. The absolute velocity of the 
blade V is given in inertial coordinates by equation (9.4). 
The relative velocity in inertial coordinates is then 
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This is transformed by equation (9.1) to give the relative 
velocity in blade coordinates 
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(10.7) 


Elements in the transfbrmation matrix and in the 
product of the transformation matrix and its time derivative 
must be kept to adequate order to produce quadratic terms in 
equation (10.7), but only those quadratic terms which 
involve time derivatives of displacements. While only 
linear terms are required for the distributed load 


SfiSSm • 

* ^ QIMUTY 


expressions, these pertleuler quadratic terms are required 
for the work expression in Section X0.3. The matrices are 
relegated to Appendix B. The required velocity components 
are then 


♦UgjCos>fp8iny ♦ cos/^cosf qj^ ♦ cos^pSiny q^ 

♦ sin/fp q^ ♦ ecos/fpSiny 

- ecos^ cosy y„ ecos^ i (10.8) 

p ‘ P t 

- Oxcos^p ♦ 

- 0(xsin^ -e)4 - Osin/9 v - tsiny ♦ y_siny3q„ 

P t P & A 

♦ [cosy - #^siny]q^ ♦ C^j^cosy ♦ 

- [{xsin^ -eMcosy - y siny) + ft xcos^ cosy wcos^ 

p Z C P p A 

- [(xsin^ -eHsiny ♦ y cosy) ♦ xcobJ/ siny ♦ wcos^ sinyjy 

p Z t p p X 

+ [xcos>9 - (xsin^ -e);? - wsin/9 3y ♦ v (10.9) 

p p t p z 

Uz - UcrC«in>fpSiny + (y^ ♦ y 2 Cosy)sin^p ^^cos^pSiny] 
u^^tcos^p - (yj^siny - y^cosy ♦ >»^)sin/fp3 
Osin^ V - (sin^ cosy - y.sin^ siny cos^ cosy]q„ 

P P ^ P w p A 

- [^in^pSiny ♦ (y^^ y^cosy )sin>fp ■»• ^^cos^pSinylq^ 

+ [cos/p - (yjjSiny - y^cosy ♦ /^.sin^plq^ 

♦ C(x-esin^ Xsiny ♦ y_cosy) - vcos/ cosy)y 

p Z p A 

- ((x-eain>fp)(cosy - y^siny) ♦ vcos/pSinyly^ 

- vsin/p y^ (x-esin/p)^^ ♦ w (10.10) 

These expressions contain all necessary terms Cor the 
distributed loads (10.3-5) and for the work expression 
( 10 . 1 ). 
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10.2 PITCH RATE OF THE BLADE 

Various angular valocitias as dotailad in Chaptar 9 
contributs to the pitch rata is about tha Y axir>r 

is about tha X axis following tha rotation, and 

0 ara about tha Z axis, and v* is about tha z axis, all in 
tha righthand sansa. Finally, and w* ara about lha y 

axis In tha opposita sansa. Equation (9.1) is usad to 
transform each of tha angular valocitias as raquirad. Tha 
total angular valocity of a point on tha blada alastic axis 
in tha daformad coordinata systam is 
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( 10 . 11 ) 


Linaar tartns ara adaquate in tha axpansion of tha angular 
valocity. Tha only componant raquirad is 

- Osin^ * Ocos/J [fi ♦ sin^. v’ ♦ cos^, w’ 1 

f p t b a 

♦ cosfi tcosy ♦ siny y ] ♦ sin^ f (10.12) 

P A 1 P “ 


10.3 VK)RX OF EXTERNAL FORCES 

Tha variation of work dona by tha aarodynamic forcas 
can ba axprassad both in tarms of blada daflactions, and in 
tarms of tha ganaralizad coordinatas. Z.i tarms of 
variations of tha blada alastic axis daflactions it is 
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MOMM. 

or mu 



»» ‘Px**.. * Py*».. * Pt'*..*** 

•pan 


where x , y , and z are deflect lone of the elaetic axil 
•a' ‘'•a' ta 

in the respective directions. In teras of variations of the 
generalized coordinates this is 


^ ®X^^X * Qy^^Y * 

♦ Qt»/, ♦ P,»q^ ♦ (10.14) 

where the P and 0 are the generalized loads of equation 
(9.14). Recall that v ■ and w ■ y^q^,? if wore bending 

modes er^ required, more generalized loads would be used at 
this p)int. 

Now, if the two sets of coordinate variations can be 
related, the generalized loads can be calculated. The 
variations are related through the velocities as in 
Chapter 4. Define the operator 

A ■ (10.15) 

Then the relationship between the two sets is given by 

(10.16) 

Equation (10.15) is tppli.d to t)>. v.loeity coaponsnts 
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(10*8-10), and the result is substituted into the first work 
expression (10.13). Then, coefficients of each variation in 
the two work expressions are separately equated. The 
generalized loads which result are 

■ Itp^^oMpCOSf - pyCsiny *► ^jcosf] - p^tsin^pCOSf ♦ 

- ^^sin/pSiny ^^cos^pCOSy 1 1dx (10.17) 

f 

• JCp^cos^pSiny * p^tcos^ - “ Pj[sin^p«iny ♦ 

sin^p(^jj ♦ f^cosy) *► ^^cos^pSinyl )dx (10.18) 

P,Ic®Mp 

- sin^p(^j^siny - ^^cosy *3^)])dx (10.19) 

Qv ■ fcp ecos^ sine - p [(x-esin^ )(cosf - ^-siny) •*• 

A ' y p y p ^ 

* ;J^xcos>J cosy wcos^pCOSy] - p^ [vcos^pCOSy ♦ 

- (x-esin^ Hsiny ♦ y cosy)))dx (10.20) 

p z 


Qy ■ J{-p^ecos/JpCOSy - [ (xsin^^-e) (siny ♦ y^cosy) * 

* /} xcosfi siny ♦ wcor^ siny] - p [vcosfi siny 
t p p z p 



(x-esin^ Hcosy - y siny)]}dx 
p z 

(10.21) 

°z ■ 

f{p [xcos^ - ^ (xsin^ -e) ♦ 

^ y p t p 

“ wain^ 3 " P vsln^ )dx 

p * p 

(10.22) 


){p ecosi p (x-esin/ )}dx 

' X p z p 

(10.23) 

p - 

V 


(10.24) 

P ■ 

w 


(10.25) 
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A nott if in order ft this point. Alaost ell of the 
('.xprttfiont in thii chapter are 9 iven in terns of v and e, 
in order to facilitate their expansion to include more 
bending nod^s. The last two equations are given in terns of 
7 ^ and and give generalized loads for these particular 
no^es; entirely analogous expressions could be written for 
any additional nodes. Then, each occurence of v or w would 
be expanded in the appropriate modes. Here, a simple one 
mode model is used for each. 


10.4 GENERALI ZEO LOADS 

To recap, at this point the aerodynamic loads 
(10.17-25) are expressed in terms of the distributed loads 
(10.3-S). The distributed loads in turn are expressed in 
terms of the velocity components (10.8-10) and the pitch 
rate (10.12). This substitution process is too complex to 
be reported in detail here. Only an outline of the algebra 
and the results are given. 

First, a simple form is assumed for the inflow and 
crossflow as follows 

'*ln * ^ rOxCOS/pCOSf 

u^^ - fOL (10.26) 

where X is the average inflow ratio, r is a linear wind 


shear coefficient , and ^ is the crossflow ratio. The use of 
node shapes and y^, and the definitions of X, t, and p 
facilitate the integration of the generalized loads. 

For a yaw alignnent angle ^ off the wind, Boa»nttai 
theory gives 

X ■ (^Vcost - 

+ [(^Vcost - j^ira) + ^2'*^ (10.27) 

and 

P - Vsin* (10.28) 

where V > (at hub height), 0 > is the solidity (rotor 

planform area/ disk area), and 6 is the reference pitch 

O 

setting. If the reference chord length is c^, then the Lock 

number may be defined as y ■ />ac L^/I , 

o b 

Chopra has gi^*en a simple formula for the wind shear in 
terms of a power law exponent p [35]. Rearranged, this is 

r ■ p(X ♦ l<ra)L/H (10.29) 

8 

where H is the hub height and p is between 0.15 and 0.4 
depending on the terrain. 

Unlike the derivation of the equations of motion, this 
develo|xnent of the aerodynamic loads requires an ordering 
scheme. Reasonable orders of magnitude for HAHTs are 


ftte- 


Tv' Tw' r/® 

X, 

c/Lf c/Lf f 


0 (.®) 

0 (.'’) 

0 (.‘) 


where < is of the order of q^A- In any coefficient, only 
the largest terns and those one half order snaller are 
retained. Coefficients of like harnonics are also ccmpared, 
and those nore than one half order snaller are discarded. 
Also, all r*, tf, and terms are ignored, effectively 
eliminating higher harnonics of the airloads. 

The contribution to the aerodynanic loads made by blade 
one is calculated by setting and q^ ■ ^i* Then, 
the contribution to the aerodynanic loads made by blade two 
is calculated by setting q^ » q^^ » % “ <J^, 2 » 

^ ■ y ir. The contributions from both blades are added to 
give the eleven generalized loads. Here again, the 
multiblade coordinates (9.15,16) for symmetric and 
antisymmetric elastic rotor modes are used. 

Finally, the generalized loads are substituted into the 
equations of motion (9.18). Many terms in the loads involve 
generalized coordinates or their time derivatives. These 
are subtracted from both sides of the equations and thus 
augment the C and K matrices. The remainder are added to 
the forcing vector Q. 
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These eerodynamic coefficients are given balov with the 
degrees of freedom in the same order as in equation (9.18). 
In these coefficients, various aerodynamic integrals are 
used, and D^. These are defined in A^>endix C. To aid 
in understanding the aerodynamic coefficients, they are 
arranged so that they are all equal to one for an idealized 
blade: flat, untwisted, uniform, and with • xA* 

^11 * OI^7{|LgX8^tl“COs2y3-|(2LyX+L^^^)ApSin2ylA^ 

-►f(2L_X-*-L )A COS2y)A^ 

4 7 6 o p ^ 

C,3 • 01j_,{i(Ljr-Lg,»^).in2r)A^ 

■ ™KrU<3L,X-L.» )U-cos2,l-iL,4 ,in2,)A 
C,. • tl^cos2f3-[(3L.X-L,F )sin2y)A 

15 b o3p '^6 5 4o ' 

[l-^cos2y3-|L,r8 sin2y)A 

lo bbjp dhO 

‘^17 ■ O'brt-JOLjX-L.JJsinr-iLj/peo,,)^ 

C, .■ 01 r(-(4L r-3L )sin2y]/L^ 

1.10 12 22 ZS'^ o ^ 

C -01 7(i(3L XL 8 )siny-iL,,^ cosfl/L^ 

1.11 b^ 3 24 23 o 3 22'^p ^ ' 

Cji ■ 0I^r{i(6L^X-3L^8^)^pCl-^COs2y34l'8^^o*^"^^*/^^ 

C 22 • OI^r{iLgX8^[l^cos2y34<2LyX^L^8^)^pSin2r)A^ 

S3 " 01^,y{-j(I*5r-Lg^F^)[l+COS2y33A^ 

[l-cos2y3-J(3L.X-L,8 )sin2y}/L 
Ss * OIj^7{i(3LjX“L^^^)[l^co82y3*jL3^p8in2f}/L 
^^26 ■ nJbr(|L^rJ_^U*co,2,]-lL3„^,in2,)/L 


Cj7 • OI^7{-3tj/^*lnf-j(3LjX-L^»^)eOif)/t 
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=38 • 
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t> 0 1 
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0Ib,ti^3X»^)/L 
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Chapter 11. HARNOHIC BALANCE SOLUTION OF 

EQUATIONS WITH PERIODIC COEFFICIENTS 

The equations of motion and aerodynamic loads derived 
in the preceding chapters form a second order system of 
ordinary linear differential equations vith periodic 
coefficients. A solution may be pursued using any of a 
number of techniques presented in the literature. Several 
of these are introduced briefly, followed by a detailed 
development of a general harmonic balance method useful for 
stability, steady*state response, and transient response 
calculations. 

Perhaps the most straightforward approach is direct 
numerical integration of the equations beginning at some 
chosen initial conditions. Stability is determined by 
inspecting the result for growth or decay. In a stable 
case, the steady-state response is found if the calculation 
is carried far enough, or if the correct initial conditions 
happen to be chosen. This method ignores Floquet-Liapunov 
Theory C53, 54, 55] which reverals the mathematical form the 
solution must take. Floquet theory has spawned three 
general families of solution techniques! 


1) Perturbation methods. 

2} Calculation of the Floquet transition 
matrix by numerical integration. 

3) Harmonic balance methods. 

The perturbation method was first developed by Hsu 
[56]. These methods are limited to cases where the 
periodicity of coefficients can be expressed in terms of a 
small parameter. Inspection shows that this is not the case 
for the equations in question. 

Calculation of the Floquet transition matrix can bi 
accomplished by a number of numerical schemes [57, 58, 59, 
60]. In general, integration proceeds over only one period 
in these methods. They come highly recommended for systems 
with many degrees of freedom. 

Hill's method of infinite determinants is a classic 
harmonic balance method [55, 61]. Bolotin applied this 
method to problems of mechanical stability, but without the 
aid of the digital computer he sought only limited 
approximate solutions to the unwieldy determinants [55]. 
Recently, Takahashi has updated the harmonic balance method 
for the stability problem [62]. The method to be developed 
in this chapter is closely related to Takahashi *s method and 
to a similar method used by Sheu [44]. Peters and Ormiston 
have derived a general harmonic balance operator for the 
steady-state response problem [63]. See also reference 
[64]. 
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11.1 THE HARMONIC BALANCE TRANSFORM 

The equations of notion and aerodynamic forces are in 
the general form 


CM(r)]{q) ♦ iCif)nq) * [K(f)Hq) - (Q(f)J (U.l) 


where new y ■ Ot and (*) ■ d/dy. The periodic coefficient 
matrices may be written as 


N 

M{y) ■ M ♦ Z (M sin ny * M cos ny) 

° n-1 ®n ^n 

N 

C(y) ■ C ♦ Z (C sin ny ♦ C_ cos ny) 

° n-1 ®n S 

N 

K(y) - K + Z (K sin ny K cos ny) 

® n-1 ®n ^n 

N 

Q(y) - Q Z (Q sin ny Q cos ny) 

° n-1 ®n *^n 


Floquet theory gives the form of the solution as [55] 


q - exp(py) ■** Z [a sin my b cos my]} (11.6) 

o m— in n 


Here, the vectors b , a and b are independent of time. 

Q B B 

A more general solution form is [see 44] 


1 



( 11 . 2 ) 

(11.3) 

(11.4) 

(11.5) 
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q • 4 q (f) ♦ Z Cq (f)sin q (f)cos mf] (11.7) 

o • s c 

m* 1 m a 

Now, as indicatad, the vectors q , q and q are 

os c 

a a 

functions of time so that 

q * ijq Z [(q ~mq )sin nf 
° m*l ®m *^m 

'** (q ‘*'mq )cos mf] (11.8) 

^m ®m 

and 

q * ^q Z [Cq -2mq -m q )sin mf 
° m»l ®m ^m ®m 

-t* (q '»’2mq *m q )cos mf] (11.9) 
®m ^m 

Equations (11. 2-5) and (11.7-9) are substituted into 
equation (11.1) and simplified using the following 

trigonometric identities. 

sin mf sin nf ■ J 5 tcos(m-n)f - cos(m+n)f] 

sin mf cos nf ■ i 5 (sin(m-n)f * sin(m+n)f] (11.10) 

cos mf cos nf ■ )s[cos(m-n)f ♦ cos(m‘*‘n)f] 

The resulting double series equation is rearranged to expose 

the sum over harmonics, and the series are truncated at some 
harmonic, P. Because this equation in harmonic series must 
be satisfied for all time f, the coefficients of each 
harmonic must balance independently. Thus, separate 
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tquatlons can ba written for conitant terns (seroth 
harnonlc), for sin f and cos f terns (first harnonic), and 
so on» up to the Pth harmonic. Generally, it will not be 
possible to balance some terms of harmonic greater than P 
which occur in the s\mis; such terns are discarded. 

This process transforms the periodic-coefficient system 
(11.1) into an approximate constant-coefficient system which 
is 2P 1 times larger. A new vector of coordinates is 
defined by stacking the harmonic coefficients, 


r ^ 



Then, the transformed equations are given in the form 

(MHq} ♦ ICHq) (RH5) • (51 (11.12) 

Bach of the barred matrices is a matrix of smaller matrices. 
For periodic coefficients up to* second harmonic (N ■ 2) and 
truncation at P *■ 3, these constant-coefficient barred 
matrices are given on the following pages. A pattern 
emerges which may be used to extend these rv<>trices for 
P > 3. 
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These equations a^ jf^^^vidCMusIng standard techniques 
for constant coefficient equations, and the periodic 
solution is reconstructed using equation (11.7)* The 
equivalence of the two forms (11.6 and 11.7) may be seen by 
realizing that the solution must be in the form 


{q} • {a}exp(pf) 


(11.17) 


where a is a vector of amplitudes. Thus, the stability of 
the periodic-coefficient system is approximately determined 
by examining the stability of the transformed equations. 
For a stable system, the steady-state response is simply 

Cq} - IR)‘45) (11.18) 


The steady-state periodic response is then determined from 
equation (11.7). It should be noted that it may be possible 
to calculate a steady response for unstable cases as well. 
The steady portion of the present method is analogous to the 
method of reference [63]. 
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11.2 INITIAL CONDITIONS AND TRANSIENT RESPONSE 

Application wf tht harmonic balanca mathod to tha 

transiant rasponsa problam is somawhat mora complicatad. 

Tha initial conditions on q and q ara not sufficiant to 

# 

datarmine initial conditions on q and ^ bacausa thara ara 
2P '•* 1 times as many. Howavar, assuming that thay can be 
established, standard techniques again apply and the 
response can be reconstructed using equation (11.7). 

To review, the constant-coefficient equations can be 
recast in state vector form as 


where 



{f} - a]{7i) - m 


(11.19) 


CA3 


0 

1 




{Rl 



The eigenvalues of the system p^ , and tha corresponding 
eigenvectors are easily computed using standard 
eigenvalue routines. Tha general solution of equation 
(11.19) is a superposition of these solutions, 


(x) - [Vj . . . ]{Cjaxp(pjf )} (11.20) 

where the c^ are arbitrary constants which may be 
determined from tha initial conditions x(0), 


(c^l 


Iv, 


.V, ...)"Mx(0)) 


( 11 . 21 ) 
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Note that the second half of each el9envector, which 
corresponds to the velocities, is no lon9er required once 
the initial conditions have been applied. The ei9envectors 
are partitioned in the same manner as q, that is 


V 

o 




(v) 



( 11 . 22 ) 






Eigenvalues p^ occur as either complex conjugate root 
pairs or as real roots, and the corresponding eigenvectors 
likewise. Since the initial conditions are real numbers, it 
can be shown that the constants c^ also follow the same 
pattern. As a result, it is convenient to combine the 
contributions from conjugate pairs when reconstructing the 
periodic response. Consider a set of conjugate pairs: 
Eigenvalue a * jLv a “ Iy 

Eigenvector u Xw U - Xw 

Constant a Xb a - Xb 


The combined contribution of 
to the response is 


such a generic conjugate pair 


21 jOA« ^ 

VTliAWP 




oTSS^ '**« • 

«Mury 


Aq ■ txp(«f) {{•U^-bl#^)CO» rf - (bU^4AI»^)tin rf 
p 

■♦• X t(-«u. ♦bw. -bu -«w. )tin 
n-1 *n *n ®n 

♦(au. -bw. -bu^ -aw^ )iin (r^n)f 
■n n n n 

♦(au "bw^ -bu *aw. )eot (r-n)» 

S S *n *n 

♦(au -bw ♦bu ♦aw^ }cot i**n)f]) (11.23) 
®n ®n *n *n 

Hera, the imaginary part of the exponential in the 

folution hai been expanded and combined with the harmonica 
in the solution. The contribution of a generic real root is 
simply 

P 

Aq • a exp(af} C^u ♦ X [u sin nj^ ♦ u cos nf]) (11.24) 

° n«l *n ®n 

The question of initial conditions is answered by 

considering the implications of applying initial conditions 

to various harmonic coefficients. An initial condition 

applied to a xeroth harmonic coefficient (q and q ) 

0 o 

implies an initial displacement or disturbance. But any 

4 

non-sero initial condition applied to a first or higher 
harmonic coefficient implies an on-going periodic motion, 
which would have to satisfy t\e equations of motion. 
Indeed, such is the case for a transient response which 
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begins at one steady-state condition and equilibrates 

at a second. Except for this type of problen, the 

conditions are applied to q^ and q as follows 

o o 


4o(0) 

q^o) 


q 




4 


«1 


24<0) 

2q(0) 


ft 

•1 


again 

initial 


r'.25) 

(11.25) 


Several notes are in order at this point. In sons 
cases, the constant-coefficient system (equation 11.11) stay 
uncouple into several smaller subsystems. A rotor with two 
identical blades, as presented in the preceding chapters, 
has two such subsystems: 

1) Even harmonics of support motion and of 
sysmietric rotor modes with odd harmonicv 
of antisymmetric rotor modes. 

2) Odd harmonics of support motion and of 
symmetric rotor modes with even harmonics 
of antisymmetric rotor modes. 

When elastic modes of the blades are included in the 
equations, these two sets are possible only if the 
multiblade coordinates are used (9.1b and 9.16). 

For stability and transient response problems, it can 
be shown that the two sets become equivalent as the number 
of harmonics P goes to infinity. One form of the solution 
is included within the other, and it is tempting to drop one 
or the other of them. However, the two subsystems are not 
equivalent when the series are truncated at some finite P. 
The subsystems will generally have different orders and 
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difftrent distributions of harsionics to th« various dagraas 
of freadon. Tha systam may ba partitionad and aaeh 
subsystan studiad saparataly, but aach subtystam must ba 
studiad a id tha rasults must ba combinad. Thasa polrits ara 
illuminatad in tha application of tha mathod in tha 


following chaptars. 




Chapter 12. THE YAW- PITCH-TEETER MODEL 


A simple rotor-tover model which may be proposed for a 
MOD-2 type wind turbine involves only yawing and pitching of 
the nacelle and teetering of the rotor. In this chapter, 
the model equations of motion are extracted from the eleven 
degree of freedom expressions derived in Chapters 9 and 10, 
and tower contributions are added. The equations are 
arranged to make use of the harmonic balance method 
described in Chapter 11. Also, the response of the 
jaw-pitch-teeter model to imbalance is calculated in closed 
form for a restricted case. In the next chapter, an 
aeroelastic stability study is presented for the complete 
yaw-pitch-teeter model. 

Thus, the yaw-pitch-teeter model is developed with 
several aims. First of all, the development demonstrates 
the transformation or reduction of the six hub degrees of 
freedom to those chosen for the tower portion of a model. 
The same process would be used to extract other models, 
simple or complex, from the eleven degree of freedom parent. 
Secondly, the yaw-pitch-teeter model is used to give 
rudimentary results for aeroelastic stability and response. 
These will help explain the aeroelastic behavior of wind 
turbines with teetering rotors. 
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Sketch 12.1 The Yaw-Pitch-Teeter Model 

Sketch 12.1 gives a schematic of the proposed model. 
The elastic motions of the tower, yaw drive, and nacelle are 
represented by an equivalent pivot of the nacelle located a 
distance t downstream from the hub point. The nacelle has 







degrees of freedom in yaw and> ih’ pitch which are 

constrained by springs and k^, and by dampers c^ and Cp 
as shown. It also has moments of inertia and Ip^ 

about the respective axes. 

The rotor has only the teeter degree of freedom 
which is constrained by a spring 2k^, and a damper 2c^» It 
spins at a constant rate 0 and otherwise has all the 
attributes introduced in Chapters 9 and 10. Tip pitch 

control is a desirable feature, and it may be included 
through the aerodynamic integrals of Appendix C. 


12.1 REDUCTION OP HUB COORDINATES 

The first step in formulating the yaw-pitch-teeter 
model is to recognize the relationship between the six hub 

deflections and the tower deflections chosen, here f and 

y 

f . This relationship is then used to transform the 

p 

equations of motion and aerodynamic loads. Quadratic terms 
are required to transform the Q matrix, since this will 
produce terms which are moved to the K matrix. 
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Sketch 12.2 Hub and Tower Deflections 


Both hub and tower deflections are shown in Sketch 

12*2. These two sets are easily related by using the 

transformations of Chapter 9. The rotations are taken in 

the same order, that is « then 4 and then f ■ Ot. In 

this case # and 4 are analogous to and 

y p XI 


4^4 ( 12 . 1 ) 

f • 0 


The vector from the pivot to the hub is given in 
inertial coordinates as (see equation 9.1) 




0^ Mqea 

^ Quautv 


0 I cos# sin# 

n ’^y ’^p 

i_ . [T^][7^]{o - 


( 12 . 2 ) 


I cos# cos# 

» n y 


Note that # and # are substituted here* This vector is 
y p 

also 




Equations (12.2) and (12.3) are compared to give 


(12.3) 


q ■ 4 cos# sin# • I f 

q - -I sin^ • -t f 

Y n y n y 

q « (1 - cos# cos# ) - 0 

Z n y P 


( 12 . 4 > 


The hub deflections q f are now all expressed in 

X z 

terms of the tower deflections f and f . The approximate 

y P 

expressions given last in equations (12.4) are adequate only 
to relate the coordinates. 

The corresponding generalized loads can be related 
through the variation of the work done by them. 
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* V'p * V’k * V’v 

*•••* 

(12.5) 

'’x ■ 



'Oy * 

* "p 

(12.6) 


- 0 



These are substituted into equation (12.5) and coefficients 
of and are separately equated to give 



Q - -t P 

X n Y 

Q + -t P 

Y n X 


‘ I f P 

n y 

- 4 ^ p 

n p 


Z 

Z 


These results could perhaps be written by inspection from 
Sketch 9.2, but the preceeding method would be vala^ble for 
more complex models. 

One way to use equations (12.1) and (12.4), and 

equation (12.7) is to set up transformations between the 

deflections and between the loads. Only the first seven of 

the eleven degrees of freedom are used in this model, 

q ,.,fi . These are related to ^ , f , and by using 

X t y p t 

(12.1) and U2.4). 
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The corresponding loads are related by using (12.7). 
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(12.9) 


Equation 12.8 is substituted into the equations of motion 
and aerodynamic terms, which are then premultiplied by the 
transformation of equation (12.9). This reduces the seven 
by seven equations to three by' three. The terms and 

need not be used when M, C, and K are premultiplied, 
but they produce linear terms when Q is premultiplied. 
These are subsequently moved into K. 
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This is a conveniant place to note an addition to the 
ordering scheme. A reasonable order of magnitude for HAin's 
is 

t A o(.'‘) 

n 

Terms of differing magnitude are combined when the equations 
are reduced, and the ordering scheme must be applied to the 
aerodynamic terms as before in Section 10.1. 


12.2 EQUATIONS OF MOTION 

All that remains after reducing the hub degrees of 
freedom is to add the tower contributions to the equations 
of motion. These are very simple diagonal mass, damping, 
and stiffness terms. Without the rotor, the nacelle 
equations of motion are 


I f 

yn'^y 

I ¥ 

pn'^p 


+ C ^ ♦ k p 

y y y y 

c / k p 
p p p p 



(12.10) 


All loads on the nacelle are ignored except for its mass 

imbalance S^. It is assumed that the system is statically 

balanced so that S • 0 and sf ■ 0. 

n bn* b 

It is convenient to define total moments of inertia 
about the yaw and pitch axes which include the rotor mass 


I 

y 




( 12 . 11 ) 


Further, the following nondimensional parameters are defined 


gl i 
YT* ! 




*:f:* 


ORIOINAL PMH m 
OF POOR QUALITY 


M 

r 


b 


y 


w2 

y 

V* 

p 


V 


2 

t 


I 

n 


k /Q^I 

y y 

k /O^i 

p p 


Ji/cVk I 

y y y 

C/L 


(12.12) 


The reduced equations of motion with aerodynamic terms 

and tower contributions are made nondimensional by dividing 

through by 20^1 , and all of the preceeding definitions are 

b 

applied. The equations are in the familiar form 
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> d/dy. The coefficient matrices are g 
pages in a form compatible with Chapter 
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12.3 AN APPROXIMATE SOZ47TION FOR IMBALANCE 

This stction is s di 9 ression from the stain dtvslopnant;. 
Tha rascwnse of the ysw*pitch*teeter model to a small rotor 
imbalance is calculated in an approximate fashion. While 
the closed form solution presented is possible only for a 
restricted case, it is nonetheless c‘ some interest. In 
this section only, the rotor is assumed to have no precone, 
no teeter spring or damper, and no aerodynamic loads. The 
in vacuo equations of motion are then 


iy ♦ ^(l-cos2f) 
-tjsin2f 
siny 


->ssin2f 


Ip ♦ ^{l+cos2y) 
-cost 


2C 1 V * sin2y 
y y y 

-1 *cos2y 
2cosf 


l-cos2y 

P P 
2slny 


2f I t iin2f 
P P P ^ 



o 

siny 


y > 


0 It* 

-tosy 


p p 

P 1 

0 0 

» 

1 

• 



( 0 ) 


(12.16) 


Imagine a small imbalance mass ft located on blade one 

at I . Besides adding negligible increments to the mass and 
81 

moment of inertia terms, general imbalance terms are 
introduced which augment X and Q, These terms are quoted 
here without documentation. They can be derived from the 
kinet^: and potentiu!. energy expressions (9.7 and 9.12) by 
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2 

replacing with m, with end with mt^ and 
applying Lagrange's equations end the transformations of 
this chapter. The resulting imbalance terms are 


and 


A[K] • (sg/n L) 


0 

0 


0 

-COSf 


0 

1 


1^0 1 cosy 


(12.19) 


A{Q} 



( 12 . 20 ) 


where s=mZ L/2I. . 

These results demonstrate what occurs whenever the 
blades are dissimilar. Compare equation (12.19) with X in 
equation (12.18). The pattern of coupling between tower 
degrees of freedom and teetering discussed in Chapter 11 has 
broken down. However, the Mathieu type terms are very small 
compared to the other stiffness terms. A small amount of 
damping would eliminate any instability they may produce, 
and they are neglected henceforth. 

The periodic excitation introduced by equation (12.20) 
is also small, but may produce resonances under some 
circumstances. These imbalance terms excite the opposite 
set of harmonics from those excited by the aerodynamic 
forcing terms (12.17). 
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An approximate solution to the in vacuo equations of 
motion (12.18) with imbalance forcin$ (12.20) is possible. 
First premultiply the equations by 


[M] 


-l 


1/ly 

0 

siny/f, 


0 

i/fp_ 

cosy/f, 


-siny/fy 

cosy/fp 

1 ♦ sin^y/fy ♦ cos^y/fp 


( 12 . 21 ) 


For the restricted case of this sect ion » y^ and yp are 
completely uncoupled and the equations become simply 







y * 2C V y + »* y - -(sZ /i^)siny 

Ty y y Ty Y ^ Y 

y * 2( y y •►r^y ■ (st /I )cosy (12.22) 
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n y p y P 

- 2coSf - 2siny y^ + •'ySiny y^ - vjcosy yp 


These equations are easily solved in sequence for y^ 

and y then fi , It C *C “0 
t y p 


y ■ 'a siny 

y y 

y > a cosy 

p P • 

^ - [a -^a ] [a -a ]cos2y (12.23) 

t p y P y 

where 

dy ■ S-tn/[ly ( «'y-l) ) 

^p ■ S^/[ly ( »'y"l) ] 


This response demonstrates an interesting mode of the 
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system. As the hub follows the imbalance, it describes an 
ellipse with radii a L and , The nacelle rotations 
are counteracted by the teetering in such a way that the tip 
path plane remains always the same. 

The results of this section apply only for the rotor 
which has no airloads, is not preconed, and has no teeter 
spring or damper. However, the rotor-tower system would 
exhibit similar behavior if these attributes are not too 
large. 
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Chapter 13. AEROELASTXC BEHAVIOR 

OF THE YAW-PITCH-TEETER MODEL 

Various aeroelastic instabilities of the teeterin9 
rotor-tower system are possible depending on the parameters 
of the machine and its operating condition. Several kinds 
of solutions of the yaw-pitch-teeter equations (12.13-17) 
are presented in this chapter for a range of parameters. 
First, transient response results are given for several 
unstable cases in the manner of Janetzke and Raza [50]. The 
results of the present study are compared to theirs as a 
means to verify the yaw-pitch-teeter model. Second, a 
stability study is r>resented which displays the effects of 
some key parameters: support stiffness, damping, inflow 

angle, and preconing. Finally, the steady response to wind 
shear is briefly discussed. 

In each solution, the NASA NOD-2 is used as a base 
case, and one or two parameters are changed at a time. The 
base case parameters are adapted from reference [50], and 
are given in Table 13.1. (Tables and Figure.* are placed at 
the end of this thesis for ready comparison.) 

13.1 COMPARISON OP TWO MODELS 

A more complete description of the model presented by 
Janetzke and Raza is in order here to point out similarities 
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and differences [503. The present model is physically 
similar to their model, although pitch is defined in the 
opposite sense. They used a flap bending mode for each 
blade in addition to the yaw, pitch, and teeter motions used 
here. They did not include preconing, but they did include 
a delta-three angle of the teetering axis. The present 
analysis assumes that the teetering axis is perpendicular to 
the blade axis. 

The aerodynamics of reference [50] included nonlinear 
expressions for the lift and drag coefficient, and so were 
not given in explicit form. This was consistent with the 
solution technique used there, which was a straightforward 
numerical integration of the equations of motion. The 
aerodynamic forces were apparently calculated numerically at 
each time step. 

The results presented by Janetzke and Kaza were 
transient response time histories, given an initial 
disturbance of the system. Stability was determined by 
examining these time histories for growth or decay, and the 
nature of the behavior was shown as well. The standard 
MOD-2 case was examined with and without damping; a case 
with reduced yaw stiffness, (v^ ■2) was examined without 
damping; and a case with both yaw and pitch stiffness 
reduced (v > v « 2) was examined with and without damping. 


144 


ORK^NAL PAGE W 
OF POOR QUALfTY 


The stability of these same standard cases may be 
examined by applying the techniques of Chapter 11 to the 
yaw-pitch-teeter equations. The harmonic coefficient 
matrices (12.14-16) are calculated and substituted into the 
barred harmonic balance matrices (11.6-8). Then, the 
eigenvalues of this transformed system are extracted and 
examined. Roots are presented in Appendix D for all of the 
cases of reference [50]. The conclusions are the same for 
the two analyses, except that the present analysis predicts 
that the standard MOD-2 case without damping is slightly 
unstable (« ■ .0015). Reference [50] had this case as 
neutrally stable. 

Appendix D gives the roots for the standard cases at 
three different truncations of the harmonic series: P * 1, 
P ■ 2, and P ■ 3. These results show the rapid convergence 
of the harmonic balance method for the stability problem. 
As P is increased, the real part of root converges and may 
repeat with an imaginary part which is different by an 
integer amount. The roots are also divided into the two 
sets discussed in Section 11.2. As P is increased, roots 
jump from set to set, but no new roots arise at P ■ 3. 

Transient response time histories such as those in 
reference [50] can also be calculated using the techniques 
of Chapter 11. The eigenvalues and eigenvectors of the 
transformed system are calculated, and the initial 
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conditions applied (11.21). Then, the contributions of all 
conjugate pairs (11.23) and real roots (11.24) are added to 
give the transient response. 

The two reduced stiffness cases without damping, both 

unatarble, are arguably the most interesting and are used 

here to compare the two models. Appendix E presents the 

eigenvalues p^ , the eigenvectors , and the combination 

constants c for the v * 2 case without damping. Figure 
j y 

13.1 presents the transient response results from the 
present analysis for this case. The initial conditions and 
all other parameters are the same as those used by Janetzke 
and Kaza, whose results are given in Figure 13.2. These are 
transformed to match the conventions used here, the bending 
mode responses are omitted, and the scale is changed. 

Figures 13.3 and 13.4 make the same comparison, but for 
the m m 2 case. Perhaps the most startling fact about 
these plots is that only harmonics up to the second are used 
(P - 2). 


The initial conditions used in reference [50] are 

apparently designed to excite the forward whirl flutter mode 

of the system. The initial conditions used in the results 

of Figures 13.1 and 13.3 are placed on the zeroeth harmonic, 

based upon the guidelines of Chapter 11. If similar 

whirling initial conditions are used but placed on the 

second harmonic, an slightly different picture emerges. 

Figure 13.5 shows this result for the v • v • 2 case, 
a y p 
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A number of unstable roots are given in Appendix D for 

the ¥ • V • 2 case. A root with « ■ 0.0104 is largely 
y p 

responsible for the phenomena in Figure 13.3, and a root 
with « ■ 0.0125 is culpable in Figure 13.5. This points 

out a difficulty with using transient response histories to 
examine stability, since the initial conditions must be 
carefully chosen to assure excitation of the unstable modes 
of interest. Thus, transient response time histories are 
.used here only for comparison to reference [50], and are not 
continued in the following stability study. 

13.2 AEROELASTIC STABILITY STUDY 

Both stability maps and plots of damping versus a 
parameter of interest are used in this section to give a 
rudimentary understanding of the aeroelastic stability of 
the model. Figure 13.6 shows the basic instability regions 
for different combinations of support stiffnesses. The 
locations of cases from the last section are shown. 
Generally, reducing either stiffness too much, or having 
them too close together can cause instability. This plot 
does not show the strength of the instability, and although 
the plot indicates that MOD-2 is unstable, a very small 
amount of damping suppresses the whole matched stiffness 
region. 

Note also that a portion of the boundary was calculated 
with both P ■ 2 and P ■ 3. This is another check of the 
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convergence of the harmonic balance method. The periodic 

coefficients contain harmonics up to second, and it appears 

that P « 2 is adequate for these calculations. 

The regions shown also must not be regarded necessarily 

as being one particular mode. Rather, many different roots 

(for P » 2 there are fifteen pairs) overlap to weave this 

pattern. This fact is demonstrated by Figure 13.7 which 

plots the real part a of various roots which are active as 

the support stiffness is reduced along the line 

Here again, the branches are separated into the two sets 

discussed in Section 11.2. This plot shows the weak nature 

of the instability near the MOD-2 base case. Note also the 

very strong instability near v « v ■ 1, 

y p 

A less confused and more rewarding picture of the 

instabilities involved is produced by decreasing the yaw 

stiffness while maintaining the pitch stiffness, as shown by 

Figure 13.8. Three general instability regions are shown: 

one at matched stiffness which is weak at least for this v ; 

P 

another beginning around • 2 which involves yaw and 
teeter as previously shown in Figure 13.1; and a divergence 

near v ■ 1. Figure 13.8 will be used as a basis for 

y 

comparison in the following plots. 

A reasonably small amount of damping applied equally to 
yaw and pitch suppresses the flutter instabilities as shown 
by Figure 13.9. Several of the roots from Figure 13.8 are 
completely off scale. 
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increasing the inflotf through the rotor disk is 
destabilizing to both the natched stiffness and reduced 
stiffness instabilities, as shown by Figure 13.10. This 
plot applies to the ease of an increase in the inflow 
without a corresponding change of the pitch setting. A 
change in the inflow which is counteracted by the pitch 
controller so that the power output is constant has only a 
minor influence on the stability. 

Finally, Figure 13.11 should be compared to Figure 
13.8, which shows the effect of a practical amount of 
downwind preconing on the stability. The matched stiffness 
instability is hardly affected, while the reduced stiffness 
instability is somewhat broader and more intense. 

The divergence region is also broadened slightly by 
precone. This last observation is perhaps a little 
surprising, and it hints at the different nature of the 
freely teetering rotor. Precone would ordinarily stabilize 
the wind turbine in yaw, but this divergence involves teeter 
as well. In this regard, see reference [65]. Stiffening 
the teetering degree of freedom changes the aeroelastic 
behavior of the machine appreciably, as shown by Figure 
13.12, which should be compared to Figure 13.6. Here, a 
large teetering spring has been applied, and the divergence 
boundary is at a much lower yaw stiffness. 




To conclude this discussion of aeroelastic behavior, 
Figure 13.13 is included. The steady response to wind shear 
(r ■ .03) is indicated for two cases of support stiffness 
with damping. The yaw and pitch response for the standard 
MOO-2 casa is actually too small to be seen on this scale. 
The teeter response is the same for both the standard case 

and the v > v ■ 2 case. 

y p 
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CONCLUSIONS 


As outlined in the introduction, this thesis has 
addressed four specific recommendations made in reference 
[14]* These same recommendations serve as an outline for 
these thoughts and are therefore repeated here. 

1) Develop simpler models to investigate the 
main origins of aeroelastic phenomena (1). 

2) Examine aeroelastic and mechanical 
instabilities more closely, especially for 
the proposed more flexible systems (3). 

3) Study teetering effects and propellor 
whirl type instabilities (4). 

4) Look in detail at generator drive train 
interaction with other system components (7). 

Two simple aeroelastic models have been developed in 
accordance with 1) above. The first, a simple equivalent 
hinge model of an isolated rotor blade, was intended to meet 
some of the requirements of 2). The second, a simple 
rotor-tower model with a teetering rotor, was intended to 
fulfill some of the requirements of 3). In the process, a 
framework was established which facilitates the develoi^ent 
of other simplified models, including possibly a model to 
address 4) . 
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In Part I, an equivalent hinge model was proposed at a 
kind of typical section approach for the isolated HAWT 
blade. The blade was assumed to be rigid, but hinged near 
the root and constrained by springs. Three degrees of 
freedom were used: flapping out of the plane of rotation; 
lagging in the plane of rotation; and torsion about the 
pitch control axis. The model derived includes offsets of 
the center of gravity and the aerodynaaiic cent radial 
offset of the equivalent hinges from the hub, and both 
precone and droop angles of the blade. 

First, the complete nonlinear equations of motion were 
derived, then they were linearized in perturbations about a 
steAdy-state blade position. The linearized equations were 
further simplified by applying an assumed ordering scheme 
reasonable for HAWTs. In particular, the ordering scheme 
allowed moderately large out of plane deflections, pitch 
setting, and inflow angle. 

Quasisteady aerodynamic loads were derived in a similar 
process. Nonlinear aerodynamic loads were written assuming 
only that the angle of attack wes small. These too were 
linearized and simplified as before. Finally, the nonlinear 
steady-state equations were written. The ordering scheme 
was also applied to these equations, which reduce to 
quadratic for flap deflection and linear for lag. The 
steady-state torsion is assumed to be prescribed by the 
power setting. 
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Th« resulting tysttm of Ilntarlxtd tquationt ond thoir 
•ssoeiattd sttxdy-statt aquations wtra iaplsnsntsd on a 
digital computar. An axtansiva paramatar variation was 
conduetad using tha NASA MOD-O KAWT as a basa easa. Tha 
stability was also calculatad for savaral two dagraa of 
fraadon subnodels axtraetad from tha thraa dagraa of fraadom 
system. These two degree of freedom models ware acceptable 
for predicting flap'torsion flutter and flap-lag 
instability, but stiff inplane instabilities ware only 
predicted by the thraa dagraa of fraadom modal, finally, 
the model showed good agreement with results from a modal 
model previously derived by the author. 

In Part XI , several tools ware first de^’alopad which 
have broad application to rotor-tower HANT 
aeroalasticity . An eleven degree of freedom linear modal of 
a teetering rotor on a flexible support was derived using 
six hub degrees of freedom: three Cartesian deflections and 
three Euler rotations. Tha use of hub motions allows the 
rotor model to be adapted to any tower modal chosen, from 
simple to complex. 

Besides teataring, tha robor blades ware modeled using 
symmetric and antisymmetric banding modes in both flap and 
lag banding; inaxtensional banding was assimad. Torsion of 
tha blades and all torsional moments ware naglactad because 
of the high torsional stiffness which characterizes most 
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HAWT blades. The rotor model derived includes preconing of 
the blades as well as a small undersling of the rotor, and 
the blades were assumed to be twisted and tapered. 

The aerodynamics were again assumed to be quasisteady. 
Whereas the in vacuo equations of motion were derived 
without an assumed ordering scheme, the aerodynamic loads 
could not have been derived conveniently without one. 
Unlike some other analyses, the aerodynamic terms were 
derived in explicit form as coefficients of the equations. 

Next, a harmonic balance method was developed to solve 
the equations of motion, which form a second order system 
with periodic coefficients. The method as outlined is 
useful for stability, transient response, and steady-state 
response calculations. The form of the harmonic balance 
method employed is quite convenient for small systems of 
equations, and the convergence has proven to be rapid. 
Thus, it would be useful as well for other problems with 
similar equations. 

The model of a teetering rotor on a flexible support 
was intended to serve as a parent from which a number of 
simpler models could be extracted. This was demonstrated by 
developing a simple rotor-tower model which includes only 
nacelle yawing, nacelle pitching, and rotor teetering. The 
resulting equations of motion were solved using the harmonic 
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balance technique. Solutions were implemented with several 
computer programs for stability, transient response, and 
steady state response. 

Transient response time histories were calculated for 
several cases taken from a study published by Janetzke and 
Raza [50]. Their model included flap bending modes of the 
blades, but did not include preconing, and they solved the 
equations of motion by direct numerical integration using 
finite time steps. The comparison between their model and 
the simple rotor-tower model developed is very good, even 
though the harmonic balance was truncated at the second 
harmonic of the rotation spee for the present analysis. 

The simple yaw-pi tch-t,eter model was used to examine 
the effect of key parameters on the whirl stability and 
divergence of a teetering HAWT. These included support 
stiffness, support damping, inflow increase due to a gust, 
and preconing. Finally, some implications of using a 
teetering rotor were investigated by including a teeter 
spring. 

In s’uiunary. the thesis contributes to both the modeling 
methodology and the understanding of aeroelastic behavior of 
wind turbines. The general rotor model shows promise for 
the development of other models which can contribute to 
understanding new aeroelastic problems which may arise. 
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Appendix A. MODAL MODEL SUMMARY 

The modal equations of reference [20] are reviewed here 
with some nomenclature changed to match that of Part I. 
Coordinate systems are the same except that, because the 
blade is twisted an angle 0.{x), so are the principal axes 9 

D 

and C* The perturbation equations are 
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+ [K] I q 
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V. 

Only the inertia and stiffness terms in the coefficient 
matrices of equation (A.l) are quoted below. Here, is 

the torsional stiffness of the control system, p, is the 
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blade material density, E is the Young's modulus, and G is 
the shear modulus. 
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Apptndlx B. 


MATRICES FOR PART XX 


This Appendix contain! tht transformation matrix (9.1), 
its time derivative matrix, and matrices which are products 
of either or both. All are expressed to adequate order for 
the expansion of the kinetic energy (9.5) and the blade 
relative velocity (10.7). The latter task is generally more 
demanding. A zero indicates that no terms significant to 
this study exist in that position of the matrix. 
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Appendix C. AERODYNAMIC INTEGRALS 

This appendix contains the integrals used in the 
aerodynamic terms of Chapter 10. Ordinarily, both the 

reference pitch setting 9^ and the reference chord c^ 
would be chosen at x/L « 0.75. 
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Appendix D. EIGENVALUES FOR MOO' 2 CASES 


This appendix presents the ei 9 envaXues calculated for 
the five MOD'2 standard cases of reference [50]. 
Eigenvalues from the first set discussed in Chapter 11 are 
given in the first column, and those from the second set are 
given in the second. The convergence of the harmonic 
balance method is demonstrated for each of these cases by 
giving the roots for truncation of the harmonic series at 
P » 1, P « 2, and P « 3. 
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Appendix E. EIGENVECTORS FOR A MOD-2 CASE 


This appendix contains tha aigtnvaluas p^ , eigenvectors 
, and combination constants c^ calculated for the ^ 

case of Figure 13.1, Only one half of each conjugate set is 
given. The eigenvalue is given first, followed by the 
associated eigenvector which is broken into harmonic parts. 
Finally, che combination constant is given at the bottom of 
the column. These constants were calculated wi;h initial 
conditions placed on the zeroth harmonic coefficient as 
follows (in radians) 


4 • 0.068 

? 

» -.136 

All other initial conditions were zero. 
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Table 6.1 STANDASO MOD-0 BLAOB PARAMETERS 
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Figure o.lO vs Stability Boundaries with = 5, 
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Figure 6.13 Root Locus Plot for Parameter ^ , Preconed Rotor 
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Locus Plot for Parameter Flat Rotor 



Figure 6.15 Root Locus Plot for Preconed Rotor 



m 

* • 

I 


tn 


I 


Figure 6.16 Root Locus Plot for Flat Rotor 
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Figure 7.1 Lag Mode Damping vs f , MOD*>0 Case of Reference [20) 
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Figure 7.3 Lag Node Damping vs ^ « MOD-0 Case of Reference [20] 
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Table 13. 

1 

STANDARD MOD' 2 PARAMETERS 


«b ■ 

5.512 

h • 1-®®® 

• 0 

^ ■ 

0.1619 

r » 0.3502 
P n 

- 0.1597 

V ■ 

y 

7.767 

Vp • 7.41 vj 

- 0 

^ • 

0 

f ■ 0 f 

- 0 

r ■ 

4.777 

X - 0.129 

- -.0524 

c ■ 

0.047 

T » 0.03 /» 

• 0 


0.9275 

6 • 0.0951 * 

0 i 

0.6312^^ 

c 

'-3 ■ 

0.9787 

e L.m 0.1059 ♦ 
0 4 

0.5529^ 

c 

^5 * 

1.0248 

Lbo - 1.0139 


Here, $ 

c 

is the 

tip pitch control setting. 
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Seven more highly damped 
modes (oc -.3) omitted. 
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Doping of various Modes vs P for V - <>- Standard ParaiMt.r. 
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Figure 13.8 Damping of Various Nodes vs p for )) * 7.41 and 
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Figure 13.9 Damping of Various Modes vs for r « / » 0.01, 12 ■ 7.41 
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